AD-A102  17%  DAYTON  UNXV  OH  RESEARCH  INST  F/O  20/% 

FAR  FIELO  CONDITIONS  IN  STEADY  FLOWS  WITH  A  FREE  STREAM  MACH  NU— ETC<U) 
JUN  81  K  6  6U0ERLEY  AFOSR-78-352% 

UNCLASSIFIED  UOR-TR-78-352%  AFWAL-TR-81-3064  NL 


ADA102174 


AFWAL-TR-  81  -3064 


FAR  FIELD  CONDITIONS  IN  STEADY  FLOWS 
WITH  A  FREE  STREAM  MACH  NUMBER  ONE 


K.  G.  GUDEKLEY 

University  of  Dayton  Research  Institute 
300  College  Park  Avenue 
Dayton,  Ohio  45469 


DTIC 

kELECTE 
k  JUL  29  1981 


June  1981 

Interim  Technical  Report  for  period  April  -  November  1980 


Approved  for  public  release;  distribution  unlimited. 


FLIGHT  DYNAMICS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


i'  'T 


NOTICE 

When  Government  drawings,  specifications ,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever ;  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS ,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


CHARLES  L.  KELLER 
Mathematician 
Applied  Mathematics  Group 
Analysis  &  Optimization  Branch 

FOR  THE  COMMANDER 


FREDERICK  A.  PICCHIONI,  LtCol, 
Analysis  &  Optimization  Br 
Structures  &  Dynamics  Division 


USAF 


RALPH  L.  KUSTER,  Jr.,  Col,  USAF 
Chief,  Structures  &  Dynamics  Division 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or 
if  the  addressee  is  no  longer  employed  by  your  organization  please  noti/vAFWAT./FTRgn. 
w-PAFB,  OH  45433  to  help  us  maintain  a  current  mailing  list”. 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations ,  contractual  obligations,  or  notice  on  a  specific  document. 


i.  kyrnsBr.; 

' 

Karl  G. 

f  Guderley 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

University  of  Dayton 
300  College  Park  Ave 
Dayton,  Ohio  45469 


Jam  MMtl  / 


>1.  CQN.TPQU.ING  OfFlCELtLAME  AND  ADDRESS -  -  —  .  - 

Flight  DynamicsLaboratory  j  j-')  / ,  ~ 

AFWAL/FIBRD  J 

Wriqht-Patterson  AFB,  QH  45433  66 


.  MONITORING  AGENCY  NAME  ft  ADDRESS^//  different  from  Controlling  Office)  IS.  SECURITY  CLASS,  (of  thie  rmport) 


I 


Unclassified 


15*.  DECL  ASSI  FI  CATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  (him  Roport) 


Approved  for  public  release,  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  the  mbitrmct  entered  In  Block  20,  II  ditto rant  from  Report) 


19.  KEY  WORDS  {Continue  on  reveree  elde  It  neceeeery  end  Identity  by  block  number) 

Transonic  Flow 
Far  Field  Conditions 
Mach  Number  1 

Hypergeometric  Functions  (Application) 
Laplace  Transform 


ABSTRACT  (Continue  on  rovoroo  elde  If  nocoeoery  m%d  Identify  by  block  number) 

The  far  field  in  a  steady  flow  with  a  free  stream  Mach  number  one 
is  governed  by  the  potential  equation  simplified  by  the  assumption 
that  the  deviations  from  a  sonic  parallel  flow  are  small.  In  a 
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form  of  a  similarity  solution.  To  formulate  far  field  conditions 
(conditions  imposed  by  far  field  to  the  computed  part  of  the  flow 
field  at  its  outer  edge)  one  considers  perturbations  to  the  , _ 
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dominant  part  of  the  distant  field.  The  particular  solutions 
of  the  linear  partial  differential  equation  so  obtained  have  also 
similarity  form.  The  far  field  conditions  express  the  requirement 
that  only  such  perturbations  are  admitted  which  die  out  more 
quickly  than  the  basic  flow  as  one  moves  toward  infinity.  The 
derivation  uses  a  Laplace  form  technique;  the  absence  of 
nonadmissible  particular  solutions  expresses  itself  by  the  absence 
of  certain  poles  of  the  Laplace  transform.  These  conditions 
are  ultimately  formulated  in  terms  of  quantities  available  in  the 
physical  plane.  One  obtains  integrals  extended  over  the  outer 
contour  of  the  computed  flow  field  which  contain  the  potential 
and  its  gradient. 
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SECTION  I 
INTRODUCTION 


In  formulating  far  field  conditions  for  steady  subsonic  flows 
one  usually  assumes  that  at  a  sufficient  distance  from  the  body  the 
flow  equations  can  be  approximated  by  the  potential  equation  linear¬ 
ized  for  the  vicinity  of  a  parallel  flow  with  the  assigned  free 
stream  Mach  number.  If  the  free  stream  Mach  number  is  increased, 
then  the  supersonic  region  which  arises  at  the  profile  extends  and 
the  boundary  at  which  one  is  justified  in  applying  far  field  con¬ 
ditions  of  this  kind  moves  to  larger  and  larger  distances.  In  the 
limiting  case  where  the  free  stream  Mach  number  is  one  this  boundary 
would  lie  at  infinity.  The  present  article  derives  far  field  con¬ 
ditions  which  are  applicable  if  the  free  stream  Mach  number  is  one. 

Far  field  conditions  cannot  be  formulated  unless  one  possesses 
analytical  solutions  by  which  the  far  field  can  be  represented.  In 
a  flow  with  the  free  stream  Mach  number  one,  the  basic  field  is  a 
parallel  flow  with  the  sonic  velocity  and  a  superimposed  expression 
which  describes  the  dominant  effect  caused  by  the  presence  of  a 
body.  This  dominant  term  satisfies  a  nonlinear  partial  differential 
equation,  but  ultimately  it  is  given  by  one  expression  in  the  form 
of  a  product  hypothesis.  Certain  parameters  which  one  encounters 
in  this  expression  (and  which  depend  upon  the  size  of  the  body) 
must  be  adjusted  as  the  flow  computations  progress.  To  this  basic 
field  perturbations  are  superimposed  which  are  sufficiently  small 
so  that  they  can  be  computed  from  linearized  equations.  Some 
perturbations  of  this  kind  increase  faster  as  one  goes  to  infinity 
than  the  basic  "dominant"  term  described  above.  The  far  field 
conditions  express  the  requirement  that  such  perturbations  are 
not  admissible. 

The  investigations  are  carried  out  in  the  physical  plane, 
even  though  the  hodograph  for  plane  flows  is  governed  by  linear 
equations.  In  the  physical  plane  two  and  three  dimensional  problems 
can  be  treated  in  nearly  the  same  manner.  An  important  aid. 


especially  in  the  three  dimensional  problem  is  the  fact,  that  the 
solutions  of  the  linearized  equations  can  be  expressed  in  closed 
form.  (Refs.  1,  2,  and  3).  The  author  found  the  approach  of 
Randall  particularly  useful. 

Using  coordinates  suggested  by  the  similarity  solution  for 
the  basic  flow  one  can  map  the  field  upstream  of  the  limiting 
characteristic  into  a  strip  of  finite  width  in  an  auxiliary  plane. 

In  this  plane  perturbations  can  be  treated  by  means  of  a  Laplace 
transform.  (This  is  the  technique  of  Reference  4.)  The  postulate 
that  perturbations  do  not  increase  faster  than  the  basic  flow  as 
one  goes  to  infinity,  in  the  physical  plane,  leads  to  the  requirement 
that  the  poles  in  the  right  half  of  the  plane  for  the  Laplace  transform 
vanish.  From  this  requirement  the  far  field  conditions  are  derived. 
The  behavior  of  the  Laplace  transform  at  some  of  the  other  poles 
is  used  to  adjust  the  expression  for  the  basic  singularity. 

Ultimately,  these  conditions  can  be  evaluated  in  the  physical  plane. 
Some  flexibility  in  their  application  is  obtained  by  means  of  Green's 
theorem.  The  report  derives  the  underlying  theory  and  the  formulae 
needed  in  the  practical  application.  The  results  hold  for  plane 
as  well  as  for  three-dimensional  flows. 
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SECTION  II 

PARTIAL  DIFFERENTIAL  EQUATIONS 


We  consider  the  flow  over  a  body  at  a  free  stream  Mach  num¬ 
ber  one.  Let  $  be  the  potential  that  describes  the  deviation  of 
such  a  flow  from  a  parallel  flow  with  the  free  stream  Mach  number 
one.  At  a  sufficiently  large  distance  from  the  body  this  deviation 
is  small  and  one  can  use  a  simplified  form  of  the  potential  equa¬ 
tion.  In  a  system  of  Cartesian  coordinates  x,  y,  z  in  which  the 
x  axis  has  the  free  stream  direction,  this  partial  differential 
equation  is 


$  + 

yy 


zz 


=  0 


(1) 


where  y  is  the  ratio  of  the  specific  heats.  In  cylindrical  coordi¬ 
nates  x,  r,  0  with  y  =  r  cos0,  and  z  =  r  sin0 ,  it  assumes  the  form 


'  <y-r>  A  (4x>  +  V  +  K  +  Vee  -  0  <2> 

Let  $o  be  the  part  of  the  potential  which  is  dominant  at 
large  distances.  In  the  three  dimensional  case  it  is  axisymmetric 
and  therefore  satisfies 


1+1  _L 
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(3) 


The  corresponding  equation  for  plane  flow  reads 

-  J-  (<b  )2  +  $  =0 

2  3x  v  o,x'  o,yy 

The  Eqns.  (3)  and  (4)  may  be  written  as  a  single  equation 


-  )2  +  -  $  +4>  =0 

2  3x  o,x  n  o,n  o,nn 


(4) 


(5) 
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where  respectively  for  plane  and  axisymmetric  flows 

a  =  0,  n  =  y  and  a  =  1,  n  =  r. 

The  solution  $  has  the  form 
o 

*Q  =  n3n"2  f(c)  (6) 

where 

C  »  (Y+1) -1/3  x  rfn  (7) 

The  constant  exponent  n  must  be  chosen  in  such  a  manner  that  cer¬ 
tain  boundary  conditions  are  satisfied.  One  arrives  at  the  famil¬ 
iar  equation 

(n2?2  -  f')  f"  +  (-5n+5-ct)n£f '  +  (3n-2)  (3n-3+a)  f  =  0 

(8) 

The  lines  £  =  -°°  and  £  =  0  are  respectively  the  negative  x  axis 

and  the  positive  y  axis.  At  the  sonic  line  one  has  <&0/X  =  0 

therefore  it  is  given  by  f'  =0.  The  coefficient  of  f",  viz. 

2  2 

(n  £  —  f * )  vanishes  for  a  value  of  C  for  which  the  generalized 

pareibola  5=(y+l)"i//3xy~n  happens  to  have  the  characteristic  direc¬ 
tion.  This  is  seen  in  the  following  manner.  The  characteristics 
of  Eq.  (5)  have  the  directions 

dx/dn  =+(y+1)1/2$„  3/2 

“  O/X 

It  follows  from  Eq.  (7),  that  for  constant  £ 


dx/dn  =  +(y+D1/'3  ncn11-1 


If 


\\ 


then 


f'  -  n2c2  -  0 


<  =  n-1(f)1'2 


One  obtains  from  Eqs.  (6)  and  (7) 


$o,x  =  (Y+l)-1/3n2n“2f  (O 


This  leads  to  the  above  result.  If  (ns  -  f')  =0  and  f"  is  bounded, 
then  the  remaining  terms  in  the  differential  equation  must  vanish. 
This  gives  the  compatibility  conditions  for  this  characteristic. 

The  line  s  =  const  for  which  this  happens  is  called  the  limiting 
characteristic.  All  characteristics  starting  at  the  body  upstream 
of  the  limiting  characteristic  end  at  the  sonic  line,  those  start¬ 
ing  downstream  do  not  impinge  on  the  subsonic  region. 

The  function  f  must  satisfy  the  following  conditions;  along 
2  2 

the  line  where  (n  t,  -  f')  =  0,  f"  must  remain  bounded,  the 
solution  $  must  be  symmetric  with  respect  to  the  negative  x  axis 
in  the  two  dimensional  case  (a=  0) ,  and  it  must  be  free  of  singu¬ 
larities  in  the  axisymmetric  case  (a=l) .  Expressions  for  f  are 
availabe  in  a  closed  form.  For  the  axisymmetric  case  we  shall  use 
the  formulation  given  in  Ref.  3.  One  finds 

n  =  4/5  for  a  =  0  ^ 

n  -  4/7  for  a  =  1 

The  functions  f(c)  form  a  one  parameter  family.  If  one  solution 
is  given  by 

f(C)  =  1(0 
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then  others  are  given  by 
f  U)  =  y3  1(Z) 


(10) 


1 

tl 

t 


with 


X  =  y-1c 


where  y  is  an  arbitrary  positive  constant. 

Perturbations  to  this  solution  are  introduced  by  setting 


$  =  <J>  +6 

o 


(11) 


and  considering  <}>  as  small  enough  so  that  terms  of  the  second 
order  can  be  disregarded.  This  leads  to  the  following  linear 
partial  differential  equation 


-  (y+1)  ,£-(<&  <b)+d>  +  (j)  =0 

1  9x  o,x  yx  Tyy  zz 


(12) 


or  in  a  form  which  combines  plane  flows  and  three-dimensional  flows 
in  cylindrical  coordinates 

-(y+D  £  <n°»0(X*x>  +  ^  <n°y  *  «  -k  <«“♦?>  *  0  <13> 

Notice  that  the  linear  operator  acting  in  this  equation  has  a 
divergence  form.  The  basic  flow  is  expressed  in  terms  of  x,  and 
n .  This  suggests  that  these  independent  variables  be  used  also 
in  the  equations  for  the  perturbations .  We  set  accordingly 


<j>  (x,n/0)  =  <J>  (n»c»0) 


(14) 


Substituting  Eqs .  (6)  and  (14)  into  Eq.  (13)  one  obtains 


(f  “  n^5^)<j>r_  +  If"  +  (-n^-n  +an)c]$  -  ani_  +  2ncn4> 

S  M  nc 

-  -  °*ee  "  0  <15> 
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In  a  first  step  toward  bringing  this  equation  into  a  divergence 
form  we  eliminate  the  mixed  second  derivative  by  setting 

♦(n»C,0)  =  ?(P,C,0)  (16) 

with 

P  =  n  h(5)  (17) 


With 


ill 

h 


nc 

f'-n2C 


one  obtains 


(18) 


(f'-nV)  (J.^  +  [f"  +  (-n2-n+an)cl^ 


f  ~-nTc2*  tp2(|>pp  +  (n+a)p<|)p]  "  a4,ee  =  0 


(19) 


Using  Eq.  (18)  one  derives  (or  verifies)  the  following  divergence 
form  of  this  equation 


~  „n+a-2.,,  2  2. 

JLrP _ (f'-n  C  ) 

3?  hn+a-l 


■*5(p»c#0)1 


n+a 


(f ’-n2c2)hn+a-1 


4>p(p*C/0)] 


-  n+a- 2 

-  °  *e  (e r t - 0)  1  -  o 


(20) 

The  solutions  of  Eq.  (19)  will  be  studied  by  means  of  a  two  sided 
Laplace  transform  for  the  p  direction.  For  this  purpose  we  set 
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r 


p  =  exp  t 

and 


(21) 


n  =  exp  t  (22) 

Then,  by  Eq.  (17) 

t  =  t  +  log  h  (23) 

One  obtains  from  Eq.  (15) 

(f  '-n2c2)  <J>  (t,5#  6)  +  If"  +  (-n2-n+an)  £]  ^(t,C,0) 

+  2n£4>r  t(t,c,0)  +  (-a+l)i  (t,5,6) 

-  4>tt(t,S,6)  -  a  ^60(t,C,8)  -  0  (24) 


and  from  Eq.  (20) 


^leXP(6t»  V^'9'1  -  A!eXp(6t> 


- — n+a-i 

(f'-n  e)hn  1  r 


-  a  JL[exp(Pt)  ~nta-l  4*0  (t, C »0)  1  =  0 
h 


with 

8  =  (n+a-1) 

Now  the  specific  form  of  the  function  f  is  introduced, 
function  f  can  be  expressed  in  the  following  parametric  form. 


(25) 

(26) 
The 
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For  a  =  0,  n  =  4/5 

1  =  -(1/2 )o_2/5  (2-o) 

f  =  (l/24)o"1/5  (o2-6a+48)  (27) 

dc/do  =  (3/10)  o~7/5(o+(4/3) ) 

The  following  expression  (for  a=l)  is  due  to  Randall. 

For  a  =  1  n  =  4/7 

1  =  -2 o“2/7  ( (1/2) -a) 

?  =  -(8/9)o1/7  (-2o2+3a+6)  (28) 

(10/7) a-9^7  (a+ (1/5) ) 

We  list  a  number  of  expression  in  terms  of  the  variable  a,  includ¬ 
ing  some,  which  will  occur  later.  Let  the  value  of  o  for  the  limit¬ 
ing  characteristic  be  denoted  by  o_ . 

Li 

For  a  =  0 


aL  -  16/3 

.  -1/2,  . -5/6 

h  =  a  '  (aL  -a) 

df/d?  =  (l/4)o1/5(o-4) 

df/d?  -  n2?2  =  -(9/100)o~4/5(a  +(4/3)  (aL 


(29) 


-  o) 


B  =  - (1/5) 

a(o)  *  a1/2(aL  -o)5/6  (30) 

b(o)  =  (25/9)  o“1/2  (4-o)  (oL  -o)''7/6 
c(o)  =  0 
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For  a  =  1 


aL  =  6/5 

h  =  a"1/2(aL  -a)"7/10 
df/dc  =  -(8/3)  03/7<-o+l) 

df/d;  -  n2^2  =  -(200/147)a"4/7 (a+(l/5) ) (aL-o) 


8  =  4/7 

a(a)  =  a(aL-a)7//5 
b (a)  =  (147/50) (1-a) (aL-a)_3/5 


c(a)  =  (3/2)a"1(a+  £) (0L-a) 


^^2/5 


The  point  £  =  corresponds  to  a  =  0.  There  one  has  for  a  =  0 


c  -  <rH)-1/3  X  y-4'5  -  -V  c'2/5  +  ... 


Hence 


a  -  (y+l)5/6  (-x/y)"5/2  y2  +  ... 


and  for  a  =  1  (33) 

a  =  (Y+l)7/6(-x/y)”7/2  y2  +  ... 

2  2 

We  observed  above  that  for  the  limiting  characteristic  (n  £  -f') 

=  0.  One  recognizes  from  the  last  of  the  Eqs.  (29)  and  (30),  that 
this  is  indeed  the  case  for  o  =  oL* 

One  uses  the  above  formula  to  replace  the  independent  vari¬ 
able  £  by  a.  The  primary  reason  is,  x>f  course,  the  fact  that  one 

then  obtains  a  simple  representation  for  f.  The  introduction  of 

1/2 

a  (more  clearly  of  o'  )  has  the  further  advantage  that  according 

1/2 

to  Eqs.  (33)  the  lines  o'  -  const  sweep  out  the  upper  half  of 
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xy  plane  from  the  negative  x  axis  to  the  limiting  characteristic 

in  a  smooth  manner  as  a  changes  from  0  to  aL.  In  contrast, c  tends 

to  negative  infinity  as  one  approaches  the  negative  x  axis.  In 
1/2 

this  regard  0  would  be  a  preferable  independent  variable,  we 
retain  a  because  of  certain  other  advantages. 

One  has  according  to  the  above  formulae 


for 

the 

negative  axis 

a 

SI 

0 

for 

the 

n  axis 

0 

= 

2; 

a 

=  0 

a 

3 

1/2; 

a 

=  1 

for 

the 

sonic  line  (f'  =0) 

a 

3 

4; 

a 

=  0 

0 

= 

1; 

a 

=  1 

for 

the 

limiting  characteristic 

a 

3 

1; 

a 

=  1 

0 

3 

°L  " 

16/3; 

a  =  0 

a  =  oL  =  6/5;  a=l 

The  values  of  a  would  change  if  one  defines  the  function  Ttlf) 
differently  (by  replacing  y  by  p  times  a  constant) .  Certain  sim¬ 
plifications  would  result  if  one  would  make  aL  =  1. 


The  coordinate  changes  described  above  amount  to  a  mapping 
from  the  xy  plane  to  the  pa  plane.  One  readily  derives  the 
following  relations: 


for  a 

=  0 

y 

=P01/2(Ol  -„)*/« 

X 

«  -(y+D1/3u  p4/5(i  - 

0/2)  (aL 

-0) 

2/3 

for  a 

=  1 

n 

=  P  01/2(o  -o)7/1° 

L 

X 

=  -(Y+1)3-/3  y  p4/7(l 

-  2a) (oL 

-0) 

2/5 

(34) 


The  inverse  transformation  (needed  in  the  application)  is  more 
complicated.  Assume  that  u  is  known,  then  from  Egs.  (7)  and  (10) 

X  «  U_1 (Y+  l)”1/3*  n“n  (35) 
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The  values  of  o  is  then  obtained  from  the  expressions  for  £  given 
in  Eqs .  (27)  and  (28).  With  these  values  one  obtains  from  Eqs. 
(17),  (29)  and  (31) 

for  a  =  0 


p  =  y  a 


-1/2 


for  a  *  1 


P 


-1/2  ,  . -7/10 

r  a  (a-a)  ' 

L 


(36) 


The  independent  variable  t  is  expressed  in  terms  of  p  by  Eq.  (21) . 

Notice  that  for  bounded  p  both  x  and  y  tend  to  zero  as  o 
approaches  the  value  a At  the  sonic  line  one  finds  for  p  =  const 
that  dx/da  =  0;  this  means  that  the  line  p  =  const  has  a  vertical 
tangent. 

Replacing  the  independent  variable  £  by  a  one  obtains  from 
Eq.  (25) 


•j^[exp(0t)  a(a)4>a(t,a,8)  ]  +  JL[exp(gt)  b(a)  <Pt  (t,a,Q)  ] 

+  ^j-[exp(8t)  c(a)4»e(t,a,0)  ]  =  0  (37) 

The  definitions  for  8,  a(a),  b(a)  and  c(o)  are  found  in  Eqs.  (30) 
and  (32). 

The  relations  to  be  derived  now,  are  helpful  in  expressing 
the  far field  conditions,  which  appear  originally  in  the  t,o,0 
system  in  terms  of  physical  coordinates.  Moreover  they  provide 
for  some  flexibility  desirable  for  practical  applications.  Let  L  be 
the  differential  operator  acting  in  Eq.  (12)  on  <(>.  Consider  a 
fixed  curve  C  in  the  x,  y,  z  space  and  two  surfaces  and  S2 


t 

‘S 


which  have  C  as  boundary.  It  is  assumed  that  the  surface 
and  its  normal  vector  can  be  brought  into  the  surface  and  its 
normal  vector  by  a  continuous  deformation.  Let  R  be  the 
volume  swept  out  by  this  deformation.  Assume  that  within  R 


L(4>)  =0  and  L(u>)  =  0 


where  w is  some  function  which  satisfies  this  differential  equation. 
Then  one  obtains  from 


M 


L(<J>  )^dxdydz  =  0 


by  carrying  out  integrations  by  part 


1lrS1“  Z1,S2 


where 


,SX  ’//’['<Y+1,#0,X<V  -*"x)SX  + 


+  (<pzu)  -  <p<oz)  ez]  .do 


and  ll,g2  is  the  corresponding  integral  for  S2.  Here,  ex,  ey  and  e, 
are  unit  vectors  in  the  respective  coordinate  directions  and  d 5  is 
the  directed  surface  element  (pertaining  either  to  or  to  S2) . 
According  to  Eq.  (39)  the  integral  1^  g  is  invariant  against  a 
deformation  of  S,  provided  that  the  boundary  C  remains  unchanged. 

For  plane  flows  the  surfaces  S  are  replaced  by  curves 
connecting  two  fixed  points.  Let  one  such  curve  be  represented 
by  x=x(p),  y=y(p)  where  p  is  a  suitable  monotonically  changing 
parameter.  Let  the  values  of  p  pertaining  to  the  end  points  of 
this  curve  be  given  by  p^  and  p2»  Then  one  has 
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(41) 


p2 

J1  mf  [-W+1)*0rx(V"^u,x)®x  +  (V  “  <>“y)®yr"dP 

Pi 


=  J  [-(y+1)  ^0fX(<Px^  -4>wx)  (dy/dp)  -  (<Jya  -  <»u>y)(dx/dp)]dp 
PI 

In  writing  down  the  corresponding  relations  in  cylindrical  6occdi~ 
hafcefe  for'  thferthree  dimensional  case,  it  suffices  if  one  uses 
axisymmetric  surfaces  S.^  and  S2.  The  vector  normal  to  these  sur¬ 
faces  is  then  perpendicular  to  the  unit  vector  in  the  8  direction. 
The  potential  for  the  basic  flow,  4>o,  is  axisymmetric,  but  the 
perturbations  need  not  be  axisymmetric.  Let  the  surfaces  or 
S2  be  given  by  x  =  x(p) ,  r  =  r(p) .  One  then  obtains  from  Eq.  (12) 

6=2ir  p2 

I2  *  f  j'  -(y+l)r  x(<l>xw  -  <|>wx) drde-r (<f>yw  -  <|>a)y) dxd0  (42) 

9=0  i 


The  expression,  Eq.  (42) ,  can,  of  course,  be  obtained  from  Eq.  (40) 
by  specialization  to  an  axisymmetric  surface. 

The  corresponding  relation  obtained  from  Eq.  (25),  again 
for  an  axisymmetric  surface,  is: 


2tt  P2 


I3  -  J  J  exp  (et)U^-  (<t>cu>  -<t>u>c)  dtd8 

6=0  £  Lh 

+  - — --_-j -  ( -<f>u>t)d;del 

(f*-nV)hn+ct  1  t  *  J 


(43) 
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where  £  and  t  are  considered  as  functions  of  the  parameter  p. 

Here  <j>  =  <J>(t,£,0)  andu=u(t ,C,0)  and  these  functions  satisfy 
Eq.  (25).  For  two  dimensional  flow  (a=  0),  the  dependence  upon 
the  third  coordinate  (here  0)  does  not  appear  and  the  integration 
over  0  is  omitted. 


Finally,  one  has  as  a  corresponding  expression  in  the  t,o,0 
system 


0=2 


V 


-ff 

0=0  p. 


exp  ( 0 1)  [a(o)  (4>aw  — <J>uia)  dtd0-b(a)  (<j>tu>-<|>a)t) dad0] 


(44) 


Eqs.  (41)  through  (44)  are,  of  course,  in  essence  the  same  except 
for  the  choice  of  the  coordinates.  All  arise  from  a  volume  integral 
of  the  divergence  of  the  same  vector  multiplied  by  w.  They  may, 
however,  differ  by  constant  factors  because  such  factors  have  been 
disregarded  in  going  from  one  form  of  the  differential  equation  to 
another  one.  One  has 

i3  =-(y+1)"1/3i1 

I3  =-y(3/10)I4  for  ct=  0  (45) 

I3  =-y( 20/21)  I4  for  a=  1 

We  shall  find  that  the  formulation  of  far  field  conditions  leads 
to  expressions  of  the  form  (44)  .  The  surface  S  represents  the 
distant  boundary  of  the  computed  flow  field,  $  is  the  perturbation 
potential,  and  oi  is  given  by  certain  particular  solutions  of  the 
equation  for  the  perturbation  potential  which  plays  the  role  of  test 
functions.  The  theory  will  be  developed  in  the  o,t  system,  while 
for  the  practical  applications  the  xyz  system  will  be  used. 

Equation  (39)  provides  flexibility  in  the  choice  of  the  distant 
boundary  surface. 
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SECTION  III 


FORMULATION  OF  A  BOUNDARY  VALUE  PROBLEM 


The  basis  for  the  future  discussions  is  Eq.  (37) .  This 
is  a  partial  differential  equation  in  the  cr,t,0  space.  In 
principle,  one  might  also  use  Eq.  (25),  but  the  introduction  of 
a  instead  of  c  simplifies  the  evaluation.  The  field  between  the 
negative  x  axis  and  the  limiting  characteristic  is  mapped  into  the 
strip  0  <  a  <  In  the  following  derivations  we  assume  that 

the  distant  limit  of  the  computed  flow  field  is  axisymmetric, 
although  the  flow  itself  need  not  be  axisymmetric.  The  map  of 
the  outer  limit  of  the  computed  flow  field  is  then  a  curve  C 
in  the  o,t  plane  which  starts  at  the  negative  x  axis  (o=  0),  and 
ends  at  the  limiting  characteristic.  Let  this  curve  be  given  by 


t  =  ts(a) . 


(For  a  field  without  axial  symmetry  the  outer  surface  of  the 

computed  flow  field  maps  into  a  cylindrical  surface  to  be  denoted 

by  S  in  the  t,  o,0  space  with  the  curve  C  as  a  directrix,,  and 

lines  parallel  to  the  0  axis  as  generators.  In  this  case,  one  has 

0  <  0  <  2tt ,  with  periodicity  conditions  at  0=  0,  and  e  =  2tt.)  The 

far  field  conditions  will  appear  as  global  relations  between  <p  and 

its  gradient  at  this  surface.  At  the  boundary  a~  0,  one  has  for 

a  =  0  (plane  flows)  the  requirement  that  the  solution  be  either 

1/2 

symmetric  or  antisymmetric  (actually  with  respect  to  o  '  ) ,  for 
the  case  as  1  (three  dimensional  flows) ,  one  has  the  requirement 
that  the  solution  be  free  of  singularities.  Details  will  appear 
in  the  course  of  the  investigation.  Boundary  conditions  at  the 
limiting  characteristic  will  be  formulated  later. 

The  far  field  maps  into  the  region  to  the  right  of  the  surface 
S,  (t=tg  (o) )  .  The  function  $  is  continued  to  the  left  by  setting 
<f>=  0. 

The  function  $  so  defined  satisfies  the  differential 
equation  (37)  in  the  whole  strip  0  <  a  <o^  except  for  the 
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surface  S  where  <J>  and  gradient  <p  have  jumps.  Let 

A4>(o,8)  =  4»(t  (a)  +  0 ,c,  9)  -  iji(tB(o)  -0,  a, 8) 

5  S 

A4>t(a/8)=4>t(ts(o)+  0/3,8)  -  <j>t(ts(a)-0,a,e)  (46) 

A4>CT ( a, 9)  =4>ct ( ts (a)  +  0,  a,  9)  -<J>a(ts(a)  -  0,  a, 8) 

There  exists,  of  course,  a  relation  between  Ad).  ,  A<{>  ,and  Ad>) 

c  o 

The  8  component  of  the  gradient  does  not  enter,  because  the 
surface  S  is  assumed  to  be  axisymmetric . 

For  large  values  of  t  the  perturbations  to  the  basic  flow 
are  allowed  to  have  at  most  the  same  order  of  magnitude  as  the 
basic  flow  itself.  The  form  of  the  basic  flow  is  given  by  Eq.  (6) . 
One  has  n  =  4/5  for  a  =  0  and  n  =  4/7  for  ot=  1.  n  is  connected 
with  t  by  Eqs.  (17)  and  (21) .  Along  a  line  o  =  constant  (which 
coincides  with  a  line  t,  =  constant)  ,  the  potential  of  the  basic 
flow  therefore  behaves  as  exp  (2t/5)  for  a  =  0,  and  exp  (-2t/7) 
for  a  =  1.  We  therefore  impose  the  conditions 

exp(-2t/s)  <j>(t,a,8)  bounded  for  t->-  +“;a  =  0 

exp(2t/7)  4>(t,a,6)  bounded  for  t-*-  +  «;a  =  1.  (47j 

For  t  <  t  (a)  one  has  <f=  0,  by  definition.  This  is  compatible 
s 

with  the  postulate 

exp(  (-2/5) -e)  t)  4>  ( t,  a,  8)  0  for  t+  -  00 ;  a=  1 

(48) 

exp(  (2/7) -e)  t)  <J>  (t,a,  8)  ->■  0  for  t+  -  »;  a=l,  e  >  0. 

The  problem  defined  by  Eqs.  (45) ,  (47) ,  and  (48)  is  always 
solvable  even  for  arbitrary  functions  A<f>,  A<J>t ,  and  A$0,  but 
then  there  is  no  guarantee  that  4>  =  0  for  t  <  t  (o)  . 

In  order  to  examine  whether  a  set  of  functions  ♦(t  (o),  a, 8) 

5 

$.(t  (a)o,8)  and  <b  ( t  (a)  ,  a,  8)  pertain  to  a  field  which  satisfies 

t  S  u  S  “ 


the  far  field  conditions,  we  set  A<p,  A#  ,  and  A$o  equal  to  these 
functions  <j>,  <f>t,  and  <J>0  and  solve  the  above  problem.  The 
function  <j>(t,  <t>t,  6)  so  obtained  will  satisfy  Eq.  (47);  these 
are  the  far  field  conditions;  but  only  if  <j>  =  0  to  the  left  of 
the  surface  t  =  t  (a)  will  the  limiting  values  of  <p,  <f>t,  and  <J>Q 
obtained  as  one  approaches  this  surface  from  the  right,  be  equal 
to  the  functions  <f>(t  ,  a, 8),  Q+lt  ,o,d)  and  A  (t  ,o,0)  which  are 
to  be  examined,  and  with  which  the  computation  started. 

Conditions  on  4>(t  .0,8),  <J>.(t_o,8)  and  <j>  (t_, a, 8)  which 

s  t  s  OS 

ensure  that  <(>  =  0  to  the  left  of  t  =  t „(o)  will  therefore  con- 

s 

stitute  the  desired  far  field  conditions.  In  the  next  section 
such  conditions  will  be  derived  by  means  of  Laplace  transform. 

In  the  case  a  =  0,  the  sonic  line  (of  the  basic  flow)  is 
given  by  o  =  4.  This  can  also  be  seen  from  Eq.  (37);  the  co¬ 
efficients  a  (a)  and  b(o)  of  <t>oo  and  4>tt  have  according  to  Eq. 
(30),  the  same  sign  for  o  <  4  (elliptic  region)  and  the  opposite 
sign  for  a  >  4  (hyperbolic  region)  .  The  direction  of  the  charac¬ 
teristics  is  determined  by  the  specific  form  of  the  functions 
a(o)  and  b(o)  .  One  obtains 


Hence, 


dt  =  +  5  (a-4)1/2 _ 

da  “  3  a1/2(16/3)-a) 

for  the  vicinity  of  a  =  =  16/3 

t  -  *  flo9  (aL~  a)  +  const. 


a*  0 


(49) 


As  a  aT  ,  t  approaches  ±  »  for  the  two  families  of 

li 

characteristics.  Incidentally,  the  integration  required  in 
Eq.  (49)  can  be  carried  out  in  closed  form 


t 


±{(10/3) log [c1/2  -  (c-4)1/2]  +  (5/3  log  [a1/2  +  2(o-4)1/2] 

(50) 


5/6  log [16/3)  -a]}  +  const. 


18 


J 


We  mentioned  before  that  because  of  Eq.  (21)  for  each 
finite  p,  which  is  equivalent  to  finite  t,  x  0,  y  0,  or  r  +0, 
as  a  0L.  Obviously,  the  outer  boundary  of  the  computed  flow 
field  given  by  the  surface  S  cannot  extend  to  the  origin  as  one 
approaches  the  limiting  characteristic.  It  will  terminate  there 
at  a  finite  value  of  y.  This  behavior  is  more  clearly  represented 
in  a  a,x  plane  (x  is  defined  in  Eq.  (22) .  Instead  of  Eq.  (37) , 
which  originates  from  Eq.  (25)  one  would  then  deal  with  the  counter 
part  of  Eq.  (24) .  One  has  for  the  characteristics  (from  Eqs.  (23) 
and  (49) 

dx/da  =  dt/da  -  dlogh/do  =  (1/3)  a"1  (16/3-cr) _1 

[±5o1/2(o-4)1/2  +  4(2  -o)]  (51) 

The  term  in  the  bracket  is  regular  for  o  =  16/3;  it  has  a  zero  of 
the  first  order  if  one  chooses  the  positive  sign.  This  cancels 
the  factor  ((16/3)  -0) -1  in  front  of  the  bracket.  For  this  family 
of  characteristics  x  assumes  a  finite  value  at  the  limiting 
characteristic . 

For  a  =  1,  one  first  carries  out  a  Fourier  decomposition 
with  respect  to  the  0  direction  and  considers  each  component 
separately.  One  then  obtains  analogous  results. 
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SECTION  IV 

LAPLACE  TRANSFORMATION 


The  problem  formulated  in  Section  III  is  now  treated  by 
a  two  sided  Laplace  transform.  The  starting  point  is  Eq.  (37). 
First  we  make  a  Fourier  decomposition  of  $(t,o,0)  with  respect 
to  9 . 

00 

<}>(t,0,e)  =  <t^m)  (t,c)  exp (im6)  (52) 

m=-°° 


Hence 


2  TT 


<f>(t,o,d)  exp(-im0)d0 


(53) 


The  Laplace  transform  of  <j>*m^(t,a)  is  defined  by 


,  (m)  , 

1 1>  (s 


,a)  =  /  ♦<"> 


(t,a)  exp(-st)dt 


(54) 


We  substitute  Eq.  (52)  into  Eq.  (37)  and  consider  only  one  Fourier 
component.  The  resulting  equation  is  multiplied  by  exp(-(s+B)t) 
and  integrated  from  t  =  -°°  to  t  =  +°°.  One  then  obtains 


t  (a)  +<=° 

/*/( 

-00  ts(a) 


~  (a(a)4>a(,n) )+ b(a)  (0(j)^m)  +  ^  ) 


(55) 


2 _ _  (m)  J 


-m  c  (a)  4> 


exp(-st)dt  =  0 


In  rewriting  this  equation  in  terms  of  ip  ^  one  must  take  into 
account  that  <p  ^  ,  <p  ^  and  are  discontinuous  along  the  curve 

t  =  tg(o).  The  individual  terms  are  considered  separately.  One 
has  the  trivial  result 
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va)  +- 


J  +  f  (-m2c  (a)  <j>  (t,a 


)  exp(-st)dt  «  m2c (o)^ ^  (s,o) )  (56) 


-<»  t  ( a) 


A  familiar  procedure  gives 


t«,(a)  +°° 


J  +  y  b(a)  (g<ji  ^  (t,a)  +  <p  (t,a))  exp(-st)dt 


(m) 


tg  (a) 


=  b  (a)  |  -  [  ( B+s)  A<t> (m)  (a)  +  A<{ >{^)(a)]  exp(-sts(a)) 


+  (gs  +  s2)  ip  ^  (s,c)  | 


(57) 


(m) 


The  notation  A<j> v  '  etc.  used  here  is  analogous  to  the  definition 

Eq.  (46)  .  For  the  evaluation  of  the  remaining  term  we  rewrite  the 

/m\ 

following  expression  in  terms  of  <j> '  ;  (t,a) 

ts  (a)  +00 

(ml  t  r 

J_ 

da 


tg(a)  +0O 

Ia(J)31,(^(S,°))  -  ^|«(o)  Js[j  .  J  *  (m)  (t,a)  exp(-st)dtj  | 


-oo  tg(a) 


ts  (a)  +°° 


9a 


9t . 


-a  (0)-g^-  A<J>  ^  (a)  exp(-stg(a)  +  a(a) 


[/  *  / 

-oo  ts(a) 


(t,a)exp(-st)dt 


^  {-a (a)  (dt  /da)A4>^  (a)exp(-st  (o)}-a(a)(dt  /da)A4  ^  (c)exp(-st  (a) 

S  S  SO  a 


9a 


+  a’  ( 


tg  (a)  +oo 

-0°  ts(a) 


ts(0)  +°o 


(t,a) exp (-st) dt 


.  +  m(aU  / 


(m) 


<^,aa*",  (t'a)exP(~st)dt 


— oo  t  (a) 

S 


] 


|-{-a(a)(dt  /da)A<fi^(a)exp(-st  (a)}-a(a)(dt  /da)A<b  ^  (a)exp(-st  (a) 
da  s  s  s  a  s 

ts  (a)  +00 


/\<J  I  -r~ 

j  To  [a(a)  <t>ain)  (t,a)]exp(-st)dt 


(58) 


ts(a) 


21 


Substituting  Eqs.  (56)  through  (58)  into  Eq.  (55),  one  obtains 

[a(a)4»  (m)  (s,a)]  +  [b(a)  (8s  +  s2)  -  m2c  (a)  ]  *(m)  (s,a) 

30  o 


=  r(m)(s,o> 


(59) 


with 

r  (m)  (s,  a)  =  b (a)  [  ( 8+s)  A<(> (m)  (a)  +  A<j>  (a)  ]  exp(-sts(a) ) 

♦  -i—  {-a(a)(dt  /da)A<f>^  (a)exp(-st  (a)}-a(a)(dt  /do)A<j>_^  (a)exp(-st  (a) 

do  S  8  8  0  8 

(60) 

Eq.  (59)  is  an  ordinary  differential  equation  for  ^  .  The 

argument  s  plays  the  role  of  a  parameter.  The  boundary  condition 

for  a  =  0  arises  from  the  properties  of  the  function  <f>  (t,a)  . 

For  a  =  0  one  obtains  the  requirement  that  \p,  considered  as  a 
1/2 

function  of  a  ,  be  either  symmetric  or  antisymmetric  with  respect 
to  a  =  0.  For  a  -  1  one  has  the  requirement  that  ij;  be  nonsingular. 

To  obtain  conditions  at  the  limiting  characteristic,  we 
consider  the  problem  in  the  t,o  plane.  One  has,  of  course, 

4> (t,o)  =  $  (t ,o) 

We  define 

+  00 

(s,o)  =  J  $^(t,a)  exp(*-ST)dx  (61) 

—00 

One  finds  with  Eqs.  (17)  ,  (21)  ,  and  (22) 

ip <m)  (s,cr)  =  hs  ij>(ro)  (s,a)  (62) 
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According  to  Eq.  (22),  r  is  finite  if  y  is  finite.  This  holds 

in  particular  at  the  intersection  of  the  surface  S  with  the  limiting 

/m\ 

characteristic.  For  finite  t  the  function  $ '  '  is  bounded  along 
the  limiting  characteristic.  The  behavior  of  <j>  postulated  in 
Eqs.  (47)  and  (48)  carries  over  to  <p  ^  and  holds  also  in  terms  of 
t.  Hence  it  follows  that  ip  is  bounded 

for  (2/5)  <  Re (6) < (2/5)  +  e,  if  a  =  0 

e  >  0  (63) 

and  for  (-2/7)  <  Re(s)<-(2/7)  +  e,  if  a  =  1 

The  behavior  of  ip  ^  is  then  found  from  Eq.  (62).  We  shall  refer 
to  the  strip  in  the  s  plane  characterized  by  Eqs.  (63)  as  the 
original  region  of  definition  of  the  Laplace  transform.  The  Laplace 
transform  in  the  remainder  of  the  s  plane  is  obtained  by  analytic 
continuation. 
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SECTION  V 

THE  HOMOGENEOUS  PART  OF  EQ  (.59) 

In  preparation  of  tha  evaluation  of  i|^m*(s,o),  we  discuss 
the  solutions  for  the  homogeneous  part  of  Eq.  (59).  We  write  this 
equation  in  the  form 

JL  (a (a)  (dG(m)/da)  +  [Xb(a)  -  m2c (o) ]G(ra)  =  0  (64) 

ao 


with 

X  =  s (s  +  6)  (65) 

Important  for  applications  is  that,  because  of  the  special  form  of 
a (a),  b (a)  and  c(o),  defined  in  Eqs.  (30)  and  (32)  this  is  a 
hypergeometric  differential  equation.  If  one  sets  in  analogy  to 
Eq.  (62) 

g<m)(a)  =  hS  G(m)(o) 


o~s/2(0l-o)"5s/6  G(o);  a  -  0 

0_s/2(0l-0)'7s/1°  G(m)(a);  a  -  1  (66) 


then  one  obtains  the  following  differential  equations 

o2((16/3)-o)  (d2g/d02)  +  ((4/3)  +  (8/3)e)0(2-o)  (dg/da) 
+  s (s-1) (0+ (4/3) ) g  =0,  a^O 


(67) 


o2(6-5o) (d2g(m)/da2)  +  12 (s+1) a (l-2o) (dg (m) /do) 

+  3(s2-m2) (l+5a)g(m)  =0,  a  =  1  (68) 

The  second  equation  is  given  in  the  article  by  Randall. 

One  readily  verifies  that 
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<t>  =  °) exP ( st) exp ( im8 )  -  g(m)  (.q)exp(sT)exp(in0) 

■  g(m*(&)y*#xp(lio0) 


(69) 


and  (because  of  Eqs.  (65)  and  (23)) 

<p  -  G(in)(,0)exp(-(s-6)  t)expCim0)  =  h"2p_6g (a)ys"6exp(im0) 

are  particular  solutions  of  Eq.  (37)  ,  and  of  Eq.  (13)  if  one  returns 
to  the  coordinates  of  the  physicaX  plane.  The  second  form  is 
singular  at  the  limiting  characteristic  if  the  first  form  is  regular, 
because  of  the  factor  h  2s  The  totality  of  the  solutions  of 
Eqs.  (67)  and  (68)  are  expressed  by  Riemann's.  P-functions, 


g  -  P 


-  0 
-s/2 


00  16/3  (=aL) 

- (s/3) + (1/3)  0 


a 


a  =  0 


( s/2)  +  (1/2)  3s 


(1/6) - ( 5s/3) 


(70) 


0 


g 


(m) 


P 


- (s/2) + (m/2) 


L -(s/2) -(m/2) 


(6s/5)+(7/10)+R 
(6s/5)+(7/10) -R 


6/5 (=oL) 

0  a 

-  ( 2/5)  -  (7s/5) 


a 


(71) 


where 

R  =  [-(3/2)m2  +  (1/4)  +  (3/50) (7s  +  2)2]1/2  (72) 


1 


Because  of  the  form  of  h  (Eqs.  (20)  and  (31)),  the  differential 
function  G  satisfies  a  hyper geometric  equation  with  the  same 
singular  points  but  different  exponents.  For  the  vicinity  of 
a  -  0,  the  exponents  of  the  differential  equation  for  G  are 
increased  by  s/2;  particular  solutions  for  G,  therefore  have  the 
form  G  *  P(o)  and  G  =  o1//2P (a)  for  a  ■  0,  and  P  *  arn/,2P(a)  and 
P  *  a  m/2P(o)  for  a  =  1.  The  solutions  for  a  =  1  which  are  free 
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of  singularities  start  with  the  term  o®^2.  Power  series  which 
start  with  the  constant  term  are  denoted  by  P(c),  used  in  a 
generic  sense . 

We  use  the  self-explanatory  notation 


g 

g 


cT*/2  § 

ysym 

-(s/2) +(1/2) 


§ 


anti 


g(m)  =  o-(s/2)+(m/2)  §(m) 


o  =  0 


a  =  l 


(73) 


G  =  (a-a)5s/6  § 
sym  L  ’sym 

r  _  ,  % 5s/6  1/2  * 

Ganti  -  loL  01  0  ^anti 

G(I”>  =  (<,L-a)7s/1°  a^2  $(m) 

Then  one  has  for  the  totality  of  the 


o  =  0 

a  =  l 

respective  solutions 


(74) 


9 


sym 


=  P 


0  °° 

0  - (5s/6) + (1/3) 

1/2  5s/2 

m 


1 

0  o/oL 

(1/6) - (5s/3) 


0 


00 


1  ' 


^anti  P 


0  - (5s/6) + (5/6) 

-1/2  (5s/2) + (1/2) 


0  o/oL 

(1/6) - (5s/3) 
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0 


00 


1 


§(m)  =  P  0  (m/2) + (7s/10) +  (7/10) +R  0 

-m  (m/2) +  (7s/10)  + (7/10)  -R  -  (2/5)  -  (7s/5) 


Now,  the  solution  with  exponent  zero  at  z  =  0  for  an  expression 


P  0 


c-a-b 


is  given  by  Gauss'  hypergeometric  series 


F(abcz)  =1+  —  z  +  a (a+1) b (b+1)  2 

F  (a  ,b,c ,  z)  1  +  cu  2  +  c(c+l)  tl  '  ” 


*  'a'u,L'4'  “  x  T  ^TT  "  T  c  (c+1)  2 1  “  ** 

Of  course,  the  notation  a,b,c  for  the  parameters  of  the  hyper geometric 
■series  (which  is  traditional)  has  nothing  to  do  with  the  functions 
a(a) ,  b(a),  c(a)defined  in  Eqs.  (30)  and  32). 

Then  one  obtains 

for  d  m  :  a  =  -(5s/6)  +  (1/3);  b  =  (5s/2);  c  =  1/2 

^  (76) 

for  ganti:  a  =  - (5s/6)  +  (5/6);  b  =  (5s/2)  +  (1/2);  c  =  3/2 

for  §(m)  a  =  (m/2)  +  (7s/10)  +(7/10)  +  R 
b  «  (m/2)  +  (7s/10)  +(7/10)  -  R 
c  =  m  +  1 

For  a  -  aT  one  of  the  exponents  is  zero,  the  other  exponent 

li 

is  given  by 


c-a-b  -  (1/6)  -  (5s/3) 


c-a-b  -  (-2/5)  -  (7s/5) 


for  a  =  0 


for  a  =  1 
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One  has  the  following  relations  between  the  hypergeometric  series 
at  z  =  0  and  z  =  1 


F(a,b,c,z)  -  r (c-a) T (c-b)  F*a'  b'  1  +  a  +  b  "  c*  1  " 


To  apply  this  relation  to  the  function  §  we  define 
^  *  F(a,  b,  c,  o/o^) 

%2  =  F(a,  b,  1  +  a+b-c,  (1  -  (a/aL))  (77) 

§3  =  (l-(a/oL)c“a-bF(c-a,  c-b,  1+c-a-b, (l-(o/oL) 

Then  one  has 


«1 


r (c) T (c-a-b)  A  .  r(c)T(a+b-c)  A 
r (c-a) r (c-b)  g2  r<a)r(b)  g3 


(78) 


We  notice  that  the  function  §2  which  is  regular  at  a  =  oL  is 
dominated  in  this  vicinity  by  §3  if  c-a-b  <  0.  This  happens  for 
Re  s  >  1/10  if  a  =  0  and  for  Re  s  >  -2/7  for  a  =  1.  These  regions 
include  the  original  regions  of  definition  of  the  Laplace  transform 
given  in  Eq.  (63).  The  same  holds  for  the  corresponding  functions 
G  (see  Eq.  (74) )  . 

In  Eqs.  (69)  we  had  used  the  functions  G  to  define  particular 
solutions  of  the  differential  equation  for  <p.  These  expressions 
can  be  used  to  obtain  some  insight  into  the  representation  of  the 
flow  field.  One  expects  that  the  function  $  will  have  the  correct 
symmetry  properties  at  the  negative  x  axis  (o  *  0)  and  that  it  will 
be  free  of  singularities  at  the  limiting  characteristic.  For  y 
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finite,  the  first  of  Eqs.  (69)  then  implies  that  at  the  limiting 
characteristic  must  be  regular  and  according  to  Eqs.  (73), 

this  implies  that  g^  and,  of  course,  also  gg  and  9ant^  must 
be  regular.  Such  regular  solutions  are  obtained  if  the  series  for 
g  terminates  after  a  finite  number  of  terms,  or  that  the  coefficient 
of  g3  in  Eq.  (78)  vanishes.  This  leads  to  the  condition  that  either 
a  or  b  is  a  negative  integer,  (provided  that  (c-a-b)  is  not  a  non¬ 
negative  integer) .  For  a  =  0,  and  gsym»  one  obtains  by  Eq.  (76) 
from 

a  =  -  ( 5s/6)  +  1/3  =  -k 

s  =  (6k/ 5  +  (2/5);  k  -  0,  1,  2... 


and  from 

b  =  5s/2  =  -k; 
s  =  - (2k/5) ,  k  -  0,  1,  2 


(80) 


for  a  =  0  and  §anti 


one 


obtains 


and 


s  =  (6k/5)  +  1;  k  =  0,  1,2... 
s  =  - ( 2k/5)  -  (1/5);  k  -  0,1,2... 


(81) 


In  neither  of  these  cases  will  c-a-b  be  an  integer.  These  values 
of  s  are  evenly  spaced.  The  spacing  for  negative  values  of  s  is 
1/3  of  that  for  positive  values  of  s. 

One  notices  that  the  square  root  R  which  occurs  in  Eq.  (71) 
for  a  =  1,  has  opposite  signs  in  a  and  b.  Therefore,  no  square 
root  appears  in  the  coefficients  of  the  hypergeometric  series  (76) . 

The  condition  for  the  termination  of  the  series  for  g^  , 
that  a  or  b  be  a  nonpositive  integer  gives 

(m/2)  +  (7s/10)  +  (7/10)  ±R  =  -k;  k  *  0,1,2... 

where  R  is  given  by  Eq.  (72).  This  is  a  quadratic  equation  for  s, 
which  ultimately  leads  to 
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k  =  0,1,2 
(82) 


s  =  (1/7) [(2k+m-l)  ±  [24k2  +  24k(m+l)  +  (6m+l) 2] 1/2 , 


A  corresponding  formula  is  found  in  the  article  by  Randall  (set 


-  1) 

• 

One 

obtains 

for 

m 

=  o. 

k 

= 

0, 

s 

=  -2/7 

and 

S  “ 

0 

for 

m 

=  o. 

k 

= 

1, 

s 

=  -6/7 

and 

s  ■ 

8/7 

for 

m 

=  o. 

k 

= 

2, 

s 

irrational 

for 

m 

=  1, 

k 

= 

0, 

s 

-  +  1 

for 

m 

=  1, 

k 

= 

1/ 

s 

=  -9/7 

and 

s  = 

13/7 

Here  it  can  happen  that  c-a-b  is  simultaneously  a  non-negative 
integral  number.  Values  can  be  found  by  a  systematic  search.  One 
has,  for 

c-a-b  =  -  (7s/5)  -  (2/5)  =  k±  kj.  -  0,  1,  2,... 
s  =  -(5^/7)  -  (2/7) 


For  these  values  of  s  one  obtains 

b  =  (m/2)  +  (1/2)  -  (kj/2)  -  (1/2)  [6 (k2  -  m2)  +  1]1/2 

For  b  to  be  a  nonpositive  integer  it  is  necessary,  but  not 

2  2 

sufficient,  that  the  radicand  6 (k.  -  m  )  +  1  be  the  square  of  an 

.  2 
integer,  say  k3.  Then 

(1/6)  (k2  -  1)  =  (1/6)  (k3  -  1)  (k3  +  1) 

must  be  an  integer.  This  means  that  k3  +  1  or  k3  -  1  must  be 
divisible  by  6. 

k3  =  6k^  +  1?  k^  *  0,  1,  2 
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Next,  one  examines  whether  there  are  integers  and  m  so  that 
k2  =  (1/6)  (k2  -  1)  +  m2 


This  is  done  by  systematically  varying  m.  Only  a  finite  number  of 
terms  need  be  examined.  If  m  is  increased,  then  the  value  of 
(obtained  from  the  last  equation  k^  increases,  but  more  slowly 
than  that  of  m.  Considering  k^  as  a  function  of  m,  not  necessarily 
an  integer,  one  has  from  the  last  equation 

dk-j/dm  =  m[  (1/6)  (k^  -  1)  +  m2] _1/2  <  1 


The  search  can  stop  after  the  difference  k^  -  m  becomes  smaller 
than  1. 


We 


k 


4 


k 


4 


k 


4 


list  some 

=  °,  k3  = 

=  !,  k3  = 


of  the  values 

+  1 ,  m=0,  k1=0  b 

m  =  1 ,  k.  =  1  b 

5,  (1/6) (k2-l)  -  4,  m  -  0,  kx  - 
7,  (1/6) (k2-l)  =  8,  m  =  1,  kx  = 
11,  (1/6) (k2-l)  =  20,  m  =  4,  k± 
13,  (1/6) (k2-l)  =  28,  m  -  6,  kj. 


=  0 
=  0 

2,  b  =  -3 

3,  b  =  -4 

=  6,  b  =  -6 
=  8,  b  =  -7 


This  shows  that  cases  are  possible  where  -b  <  c-a-b  and  where 
-b  >  c-a-b.  One  has  in  the  basic  flow  (from  Eqs .  (6)  and  (9)) 


=  y2^5  f{£)  for  a  =  0 

$0  =  r~2/^7  f(C)  for  a  =  1 


(84) 


One  will  postulate  that  none  of  the  perturbation  will  dominate 
the  basic  flow  as  y  +  ®  or  r  +  ».  According  to  Eq.  (69) ,  one  must, 
therefore,  exclude  for  a  =  0  particular  solutions  with  Re  s  >  2/5, 
and  for  a  =  1  particular  solutions  with  Re  s  >  -2/7. 
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A 


The  role  of  the  particular  solutions  for  negative  values  of 
s  is  not  obvious.  There  is,  for  instance,  the  question  whether  the 
particular  solutions  with  the  special  values  of  s  obtained  here  are 
suitable  to  represent  the  propagation  of  a  singularity  along  the 
limiting  characteristic.  Such  a  singularity  would  arise  if  it  is 
present  in  the  profile  shape.  The  discussion  of  the  problem  by 
means  of  the  Laplace  transform  which  we  have  chosen  here  will 
answer  questions  of  this  kind. 

Some  of  these  particular  solutions  have  simple  interpretations, 
which  are  important  for  the  applications. 

The  particular  solution  for  s  =  0  gives  $>  =  const,  this  is 
certainly  an  admissible  (although  uninteresting)  perturbation. 
Particular  solutions  for  a  =  0  and  s  =  1  give  4>  =  y,  this  trivial 
solution  is  excluded  because  it  alters  the  boundary  conditions  at 
infinity. 

An  expression  (for  a=  0) 

$  *  (y  -  Ay2/5(y  +  Ay)3f(!f)  (85) 

with 

X  =  (y  +  Ay)-3"(Y  +  l)-1/3(x  -  Ax)  (y  -  Ay)  ”4/5 

is  a  solution  of  the  original  problem  Eq.  (1)  with  4>  omitted,  for 

z  z 

it  is  the  expression  (10)  with  y  replaced  by  y  +  Ay  in  a  system  of 
coordinates  in  which  the  origin  is  shifted  by  Ax  in  the  x  direction 
and  by  Ay  in  the  y  direction.  Considering  Ay,  Ax,  and  Ay  as  small 
and  setting  (as  before) 

X  =  y-1?  =  y-1(Y  +  l)-1/3x  y"4/5 

one  obtains 

4>  =  Cy  -  Ay)  2//5  (.y  +  Ay)  3f  (.If) 

=  y2/5y3TCt)  +  Ay  y2  y2/5 131(c)-  ?(df/d?)] 
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-  Axy2(Y+l)‘1/3y"2/5(df/dC) 

-  Ayy3  y"3/5[(2/5)f  -  (4/5)  C  (df/dO  ] 

+  higher  order  terms 

It  follows  that,  perturbation  solutions  with  s  =  2/5  express  a 
change  of  y  (the  intensity  of  the  basic  perturbations) ,  symmetric 
particular  solutions  with  s  =  -2/5  and  antisymmetric  particular 
solutions  with  s  =  -3/5  represent,  respectively,  a  shift  of  the 
origin  in  the  x  and  y  directions.  Using  the  Eqs.  (27)  and  (29), 
one  can  express  these  functions  in  terms  of  a  and  identify  them 

with  the  functions  d _  and  d _ . .  defined  above.  One  obtains, 

for  instance,  using  the  definitions  (73) 

3fU)  -  ?(df/d?)  =  5c"1/5  =  5  "1/5gsym(°,  2/5)) 

where  the  second  argument  of  g  refers  to  the  value  of  s.  Therefore, 
ultimately 

9  =  (y~Ay)  2/5  ( y+Ay )  3f  ( zf) 

=  y2/5y3f(c)  +  Ayy2  5y2/5o_1/5§  (a,  (2/5)) 

sym 

+  Ax(Y+ir1/3u2  y‘2/501/5«syln(a,  -(2/51)  (86) 

-  Ay  u3(1/2,  y-3/5o(1/2'+(3/10>«ant(a,  -(3/5), 

+  •  •  * 

Analogous  discussions  are  carried  out  for  a  =  1.  One  starts 

with 

$  =  (y+  Ay)3{  r  +  Ar)  2/^f(t)  (87) 

with 

f  =  (y  +  Ay)_1(Y+  l)"1//3(x  -  Ax)  (r  +  Ar)  _4'/7 

and 

(r  +  Ar) 2  =  (y  -  Ay) 2  +  (z  -  Az) 2 
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setting 


?  =  M-1(Y+  1)"1/3  x  r“4/7 


Then  one  obtains  by  considering  Ay,  Ax,  Ay,  and  Az  as  small 


$  =  y3r_2/7f(?)  +  Ayy2  r‘2/7(3f  -  c(df/dO) 

-Axy2(y +  1)"1/3  r~6/7(df/d?) 

+y3r-9/7((2/7)f  +  (4/7)?(df/d?)  (Ay  cos  9 +  Az  sin  0) 
+  •  •  • 

$  =  (y  +  Ay) 3  (r  +  Ar)~2//7  f(£) 


=  y3  r“2/7  f (?)  -  Ayy215  r“2/7a1/7# ( 0)  (a  ,  -  (2/7)) 

+  Axy2(Y+D“1/3(8/3)r"6/7a3/7d(0)  (a, -(6/7)) 
-  y3(16/3)r~9/7c8/7§(1) (a, -(9/7)) 


(Ay  cos  0  +  Az  sin  0) 


(88) 


These  equations  can  be  used  to  identify  the  modifications  of  a 
basic  field  (that  is,  the  changes  Ay,  Ax,  Ay,  and  Az)  as  one 
proceeds  from  one  iteration  to  the  next  one.  Eqs.  (86)  and  (87) 
are  rewritten,  using  Eqs.  (74),  (29)  and  (31).  One  obtains 
from  Eq.  (86) 
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<J>  =  (y  -  Ay)2//5(y  +  Ay)3  f(C)  (89 

=  y2//5y3fU)  +  Ayy2  5  Gsyin(0  *  2/5)  exp  (2t/5) 

+  yxy2(Y+l)‘1/3  GgyTn(o,-(2/5))exp(-2t/5) 
-  yy  3 (1/2)  Ganti(a,-(3/5))exp(-3t/5) 


+  •  .  . 

and  from  Eq.  (88) 

$>  =  (y  +  Ay) 3 (r  +  Ar) ~2/1  f  (f)  (90) 

=  y3r"2/7f(?)  -  y2  15  G(0) (o,-2/7)exp(-2t/7) 

y2(Y+l)"1/3(8/3)G(0) (o , -6/7) exp (-6t/7) 

-  y3(16/3)  G  (a,-9/7)  exp  (-9t/7)(Ay  cos0+Az  sin0) 

The  Wronskian  of  the  functions  G^  and  G2  (for  all  cases)  is 
defined  by 

W(G1G2)  =  W(s,0)  =  (9G/3a)G2  -OG2/9o)G1  (91) 

The  following  relation  is  a  consequence  of  the  general  theory  of 
ordinary  linear  differential  equation 

a(o)W(s,a)  =  Ms)-1 

where  a (a)  is  one  of  the  coefficients  occurring  in  the  differential 
Eq.  (64)  for  the  functions  G.  Important  is  the  fact  that  this 
product  depends  only  upon  s.  The  expression  k(s)  will  appear  in 
future  formulae.  For  a  derivation,  one  writes  Eq.  (64),  once 
for  G^,  a  second  time  for  G2,  multiplies  the  first  equation  by  G2, 
and  subtracts  from  the  result  the  second  equation  multipled  by  G^. 

One  obtains 
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(92) 


which  leads  directly  to  Eq.  (90).  To  compute  the  function  k(s), 

one  needs  the  functions  and  G2  and  their  first  derivatives  with 

respect  to  a  at  one  station  a,  for  instance  at  o  =  o. .  One  has, 

•Ll 

because  of  Eqs.  (74) 


W«W'  GSym.2>  *  <°L  -  °>5S/3  "<W  «sym,2» 


W (Ganti , 1 '  Ganti , 2 


)  =  (0L  -  0>5S/3a«(<ianti>1,  «anti>2) 


w 


G<m))  =  (cL  -  a)73'5  a”  9<m)) 


(93) 


Because  of  Eq.  (78)  ,  one  has 


"<«!'  «2»  -  r'fii‘rn  '  C>’  W(«3-  «2>  «94> 


This  is  now  evaluated  for  the  vicinity  of  a  -  aL-  One  obtains 
using  Eqs.  (77) 


W(§1,  ^2)  =  -aL1  (1 


(a/aL))c  a  k  ^(c-a-b) 


r (c) T (a+b-c) 
T  (a)  r  (b) 


(a+b-c) _ _ ^  c-a-b-1  r (c) r (1+a+b-c) 

L  '°L  ~  0  r  (a)  T  (b) 


(95) 


Using  the  values  of  c-a-b  listed  before  Eq.  (77) ,  and  the  expressions 
a(a)  to  be  found  in  Eqs.  (29)  and  (30),  specialized  to  the  vicinity 
of  o  =  oL,  one  finds  from  Eq.  (92) 
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k 


sym 


(s) 


(1/3) - (5s/3)  T(a)r(b) 

l  rTcTfTT+a+b^T 


a  =  0 


<«)  =  a -(4/3)-(5s/3)  r(a)r(b) 
anti'  L  He)  r  (l+a+b-c) 


a  =  0 


k(m) 


(s) 


,-m-(7/5)-(7s/5)  T(a)r(b) 

L  r('c)  r (l+a+b-c)  ; 


a  =  1 


(95) 


where  the  pertinent  values  of  a(s),  b(s),  and  c(s)  are  lietod  in 
Eq.  (76).  In  Eqs.(79)  through  (82),  the  values  of  s  for  which 
a  or  b  are  nonpositive  integers  have  been  listed.  For  values  of  s 
in  the  vicinity  of  point  sp  where  either  a  or  b  is  a  nonpositive 
integer — k2»  one  has 


^ 

(-1)2 

+  0(1) 

i  \a/  ■ 

s  -  SP 

(da/ds)  f  (l+k2) 

r  ( —  * 

(~)k2 

■  +  0(1) 

s  -  sn 

(db/ds)  T  ( l+k_) 

(96) 

(97) 


With  these  formulae,  the  dominant  part  of  k(s)  can  be  determined 
provided  that  none  of  the  other  factors  vanishes. 

Let  sp  be  a  point  at  which,  simultaneously 

b  =  "3^2 r  k2  =  ^ ^ /  2... 
c-a-b  =  e=k^,k^=l,  2,  ... 


Let  for  neighboring  points 


b  =  -k2  +  (db/ds)  (s-Sp)  =  -k2  +  Ab 


c-a-b  =  e  =  +  Ae 


Then  one  obtains 


T(b)  _  Ae , 
F  (l+a+b-c)  “  Ab' 


k1-k2  Hk^ 
r (l+k2) 


k^  —  1,  2,  ... 
k2  =  0,  1,  2,  ... 


T(b)  =  _1_ 
f (l+a+b-c)  Ab 


rd+k2) 


0 

0,  1,  2... 


(98) 


(99) 


37 


One  obtains  from  Eg.  (78) 


_  r  (c)  T  (c-a-b)  t,, 
gl  ~  r(c-a)r(c-b)  g2 


a  =  -k2,  or  b  =  -k2 
)«2  0 ,  2 ,  ... 

c-a-b  ^  k^ 
k^  =  1 ,  2 ,  ... 


The  formula  is  applicable  also  for  c-a-b  =“k^,  k^  =  2,  3, 

b  =  “k2/  k2  =  0,  1,  2,  ...  (a  case  which  does  not  occur),  but 
a  limiting  process  is  necessary,  because  T(c-a)  =  r (-k^+b)  = 

r(-kl'k2) * 

If  at  a  point  s  =  sp,  c-a-b  =  0,  then  one  has  c-a  =  b,  i 
Therefore, 

§2(a,  Sp)  =  §3(a,  sp)  =  F(a,b, 1, 1,  l-(o/oL),  c-a-b  =  k1=0 

If,  simultaneously,  c-a-b  =  e  =  0 ,  and  b  =  -k2,  k2  *  1»2, 

then  §1#  §2,  and  §3  are  polynomials.  One  obtains 

h<°’  sp)  -  <-)k2  r(i+k2)  4{ff  92<°.sp>,  ^  ■  0 

K  2“  v  |  X  f  4  f  •  •  • 

Eq.  (102)  is  not  the  limiting  form  of  Eq.  (100) . 

Now  we  consider  cases  where  c-a-b  =  e  =  k^,k^  =  1,  2, 
and  b  =  -k2,  k2  =  0,  1,  2,...  Specializing  the  second  of  Eqs 
we  define  for  k2  <  ^ 


g2(a,  sp)  =  lim  F  (a,-k2,  l-k1+  e,l-(o/oL)) 
e-*-0 


For  e  ji  0,  and ,  therefore,  also  in  the  limit  e  Q,  the  function  F  is 
a  polynomial  which  terminates  with  the  power  (l-(a/oL))  2. 
Alternatively,  we  can  write 


g2(o,  V 


k2 

£ 


r  (a+i) 
T(a) 


r(i+k2) 
r (i+k2-iy 


r  (kj^-i)  2 

Nk1)  rd+1) 


k2  <  kl 


(100) 

•  • 

then 

-b  =  a. 

(101) 

J 

•  •  / 

(102) 

•  • 

(77)  , 


(103) 

_ J 
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IF" 


We  define,  for  k2  >  k^  as  g2,  the  same  sum,  but  terminating  with 


the  power  k^  -  1 


*1-1 


g2(°'  SP> 


r (a+i)  r(1+k2}  nki-i)  i  _  o  x i  „ 

"Z^  r (a)  '  r(i+k2-i)  r(k.  1  ru+i)  a  * 
i=0 


k2  >  kx 


(104) 


One  obtains 


r(c)T(k1) 


Si(a'  sp)  -  r  (k1-k:2)  T  (K^^+a)  g2(a'  SP)  ' 


c-a “b  f  If  2, • • . 

b  =  "k 2 #  ^2  =  0  #  1  *  2  f  • • 


k2  <  kx 


(105) 


One  obtains,  for  k^  <  k2,  from  Eq.  (77) 

kl 

§3(0,  sp)  =  (i-(o/oL))  x  F(k1-k2,  kx+a,  l+klf  (l-(a/aL))  (106) 

kl 

This  is  a  polynomial  which^begins  with  the  power  (1  -(a/o^J)  and 

terminates  with  (l-(a/cL))  l.  Here  one  obtains 

kx-k2  r(c)r(l+k2)  c-a-b  =  klf  kx  =  1,  2,  ... 

§l(o,  sp)  =  (-)  naTTTT+k^T  ^3 '  b  =  -k2,  k2  =  0,  1,  2,... 

(107) 

Notice  that  the  coefficient  of  §3  is  not  the  limit  of  the  coefficient 
of  §3  in  Eq.  (78) . 

Also,  needed  are  expressions  for  §2  in  the  vicinity  of  points 
s  =  sp  where  c-a-b  =  k^,  k^  =  0,  1,  2,  ...  The  case  k^  =  0  is 
covered  by  Eq.  (101) . 

If  b  ^  -k2,  k2  =  0,  1,  2,-..,  one  Obtains 
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kr1 


r(k,-i) 


Y-'  ,  r  (a+i)  I\(b+i)  V~1  1_ 

g2(ofs)  =  2^  (_)  r(a)  r(b)  r  (k.J  F7T+3X  (1  (a/aL}) 

i=0 


(-) 1  z 

i=k. 


s-s 


d  "P  r (a+i)  T(b+i)  r(l+Ae)  r(l-Ae)  (  1 

ds  )  Ae  r  (a)  T(b)  r(k1+Ae)  m+i-k^Ae)  J  f (i+1) 


(1 -(a/Oj) ) 


(-) 


d  1 


±  na+k.^  r  (b+kx) 


§7(o 


^  de/dis  ~Tai  tot"  rajmr+iqr  y3w  '°p 


k  =  1 ,  2 ,  3 ,  ... 
b  ^  “k ^ ,  k^  =  Of  2 , .  .  . 
c-a-b  =  e  =  k1  +  Ae 

(108) 


A3 1  functions  of  s  are  to  be  evaluated  for  s  =  sp.  Only  the 
general  structure  of  this  formula  will  be  needed.  For  the 
evaluation  of  derivatives  (not  needed  in  the  present  context)  we 
note  the  following  formulae. 


d_ 

da 


log 


r(a+i) 

Ida) 


1 

a+i-1 


f  l  8  1|  2 ,  3  f 


_d_ 

ds 


s_s 


_ P 

Ae(s) 


_d2Ae/ds2  f  c-a-b  =  k.  +  Ae 
2(dAe/ds) 


(109) 


In  cases  where »  in  addition,  b2  -  ~^2'  k2  ~  ^ 2'  ***  one  ^as 

k. 

g3(o,  sp)  =  (1-  (a/aT ) )  1  F^+a,  kj-kj,  1+^,  (l-(a/cL))  (110) 


For  k3  >  k2/  this  is  an  infinite  series,  for  kj^  <  k2,  it  is  a 
polynomial.  One  has 
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lim  §_(o,s)  =  g-(o,sD) 


s-*s_ 


+  (-) 


kl"k2 


Ab  r(a+k1)r(l+k2)r(k1-k2) 
lim  (Ae)  r  (a)  r'(k,  VrTl+k,") 


b  =  -k2  +  Ab,  c-a-b  =  e  =  k1  +  Ae,  k2  <  k^ 


(111) 


The  function  g2  is  defined  in  Eq.  (103) .  The  factor  of  is 
bounded.  Finally, 


§2  (a ,  s)  =  g2  +  Y] 

i=k. 


(-) 


kl-i 


j  /S-S 

A-( _ R 

ds\  Ae 


r(a+i) 

“rTaf" 


r (l+Ae)  T(l-Ae)  r(l+k2-Ab)  ^ 

r  (kL+Ae)  r  ( 1+i-k^-Aej  r  (l+k2-Ab-i)/ 


rTTflT  (1-(°/°Ln 


r (a+k  ) 


r(i+k_) 


+  s-s  de/ds  r  (a)  r  (k.^)  r  (l+k^k^)  r  (1+k^  ^3  ^a'sp* 


,  k2“kl  Ab  r(1+ki) 

+  (-)  a¥ 


t^L  y  mail  rll~k2)  i  a.la/a  pi 

r(k.)  La  r(k,)  ra+i-k,)  r(i-i)  10/0l;; 

i=k2+l 


kj^  s  1 1  2 ,  •  •  •  f  kj  —  0^  X f  2 •  •  • 

k2  >  kl  (112) 
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SECTION  VI 

EVALUATION  OF  THE  LAPLACE  TRANSFORM 

Now  we  return  to  Eq.  (59)  .  Let  G^,  be  a  particular  solution 
of  the  homogeneous  part  of  this  equation  (that  is  of  Eq.  (64)) 
which  satisfies  the  conditions  at  o  =  0.  The  definition  is  found 
in  Eqs.  (74 ) . 

The  boundary  condition  at  a  =  aL  is  derived  from  the  require¬ 
ment  that  in  the  original  region  of  definition  of  the  Laplace 

/_\ 

transform  ip '  1  is  bounded  (see  the  discussion  preceeding  Eq.  (63)). 

The  behavior  of  the  function  is  partially  determined 

by  the  inhomogeneous  part  r ^  defined  in  Eq.  (60)).  Assume 
momentarily  that  r ^  e  0  in  the  vicinity  of  a  -  oT .  In  this  region 
<pv  '  is  then  given  by  a  combination  of  two  linearly  independent 
functions  G,  of  G2  and  G3,  say.  According  to  Eq.  (66)  G  =  h~sg, 
therefore  ip(m^(see  Eq.(62))is  a  linear  combination  of  g2  and  g^. 

In  the  original  region  of  definition  of  the  Laplace  transform 
(given  by  the  inequalities  (63) )  ,  the  function  g3  tends  to  infinity 
at  o  =  ctl  (this  is  seen  from  Eqs.  (89)  and  (77)).  At  o  *  aL 
the  particular  solution  G3  is  therefore  not  allowed  to  be  present. 

On  this  basis  one  finds  for  the  solution  of  Eq.  (59) 

a 

ip(m)(s,a)  =  k (s)  [G1(s,a)  J  r(m)(s,a)  G2(s,c)da 

a=o 

L  (113) 

a 

-G2  (s,a)  J  r  (s^cjGj^  (s,a)dcj 
a=0 

The  conditions  at  a  =  0  and  a  =  aL  are  satisfied  by  the  choice  of 
the  lower  limits  of  the  integrals.  In  the  first  place  this 
expression  holds  in  the  region  of  definition  of  the  Laplace 
transform.  By  analytical  continuation  it  is  then  extended 
throughout  the  complex  s  plane. 
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Here  the  expression  r (Eq.  (60))  is  substituted.  This 
expression  consists  of  three  summands.  For  the  second  summand  an 
integration  by  parts  is  carried  out.  The  contributions  of  the 
terms  at  the  limits  of  the  integrals  a  -  0,  and  a  =  vanish,  the 
first  one  always,  the  second  one  for  values  of  s  in  the  original 
region  of  definition  of  the  Laplace  transform.  With  some 
rearrangement  of  the  terms  one  then  obtains 

4>(m)(s,o)  =  k(s)  J  G1(m)(s,c)K(^(m)(G2(in!aL(o)- 

G2  (m)(s,a)  K (  4>(m)  0 , a)  (114) 

where 

K(<|>(m,.G(m),a,b) 

b 

=  J  exp(-sts(a)  {<|>(in)[(e  +  s)b(a)G(m)  +  (dtg/da)a(o)  (3G(m)/3a)  ] 
a 

+  G(m)  [b(a)<t>t(m)  -  a(o)  (dtg/do  )(pg(m)]  }  da  (114a) 

Here  the  arguments  of  <j>  ^  and  its  derivatives  are  tg(5)  and  a. 

The  arguments  of  G^  and  its  derivatives  are  5  and  s.  We  set 

(t,cr,s)  =  G^  (s,a)  exp(-(B  +  s)t)  (115) 

aj(m)  (t,a,0,s)  =  aj*m)(t,a,s)  exp(-im6)  (116) 

The  function  is,  according  to  Eq.  (69)  ,  a  solution  of  the 

original  partial  differential  equation  (37)  in  the  t,  a,  0  space. 
Expressing  4>  ^  in  terms  of  $  by  Eq.  (53),  one  then  obtains 
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3* 


2ff  b 

KW(ra),GW(a(b)=-(l/2  irlf  f  expC6tl|b(3}  -  (J^ui^JdS 

9=0  a 


-  a(5)  [<j>u)im)  -  4>gu)(m)  ]dt  J  d6  (117) 

Here  it  is  assumed  that  the  path  of  integration  from  a  to  b  in 
the  5,t  plane  is  given  by  t  =  tg(5).  The  subscript  s  of  t  has 
been  omitted;  dt  stands  for  (dtg/d5)da. 

The  integral  K  has  the  form  of  Eq.  (44)  and  Eq.  (44)  is 
related  to  the  expression  (40)  ,  (41)  and  (42)  by  Eqs.  (45) . 
Accordingly,  K  can  be  evaluated  in  the  physical  space.  For  a 
closer  study  we  express  <u  ^  as  far  as  possible  in  physical 
coordinates.  According  to  Eqs.  (21)  and  (22),  one  has 

exp  t  =  yh,  a  =  0 
exp  t  =  rh,  a  =  1 

With  Eqs.  (29)  through  (32) ,  and  (74)  one  then  obtains 

V  -  y(1/5)~S  o-<1/10)+(s/2)(aL_o)(5s/3>-(l/6)§sym(o) 

“anti  -  y(V5)"S  a(2/5)^^)(aL-a)(5s/3)-(l/6)§anti(a) 

-(m)  =  y-  (4/7)  -s  0  (m/2) +2/7)  + (s/2)  (  ff)  (7s/5) + (2/5) § (m)  (<j) 

(118) 
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These  formulae  express  w  in  terms  of  y  and  a.  At  a  ?  aL,C  is  a 
regular  function  of  a,  and  £  is  expressed  hy  coordinates  in  the 
physical  plane.  It  follows  that  in  the  vicinity  of  the  curve 
o  =  oL  and  for  y  /  0  or  r  /  0(  there  is  a  one  to  one  mapping  from 
the  c, y  plane  to  the  x,  y  or  x,  r  plane.  The  singular  behavior 
of  the  functions  to  at  a  =  which  is  evident  in  Eqs.  (118) ,  there¬ 
fore  carries  over  to  the  physical  plane.  The  integration  in  the 
expression  K  (in  the  form  of  Eqs.  (40),  (41),  or  (42))  will  be 
extended  over  the  outer  boundary  surface  S  of  the  computed  flow 
field.  One  can  therefore  choose  x,  y  (or  r)  as  functions  which 
are  regular  at  a  =  The  function  <j>,  <J>^_  and  <f>G  will  be  regular, 

unless  a  singularity,  induced  by  the  shape  of  the  body,  propagates 
along  the  limiting  characteristic.  Singularities  in  the  integrands 
of  the  expressions  (40) ,  (41) ,  and  (42)  therefore  appear  solely 
because  of  singularities  in  to.  The  functions  §2  are  regular  at 
a  =  aT .  If  one  substitutes  those  into  Eqs.  (119),  then  one  finds 
that  the  strongest  singularities  in  Eqs.  (40),  (41),  and  (42)  those 
caused  by  and  wy  have  the  form  (aL~o)  ^5s/3^  “  for  a  =  0, 

and  (aL-a)-f3/5)+  ^7s/5)  fora*  1.  If,  on  the  other  hand,  the 

functions  cu  are  formed  with  (Eq.  (77)),  then  the  functions  to 

c—a— b 

are  regular  because  of  the  factor  (l-(a/aL)  in  §3*  The  same 

result  can  be  found,  from  the  expression  (117)  in  the  S,  t  plane, 
but  the  discussion  is  quite  cumbersome. 

Poles  of  the  expression  (114)  arise  for  several  reasons: 

/m\ 

a.  In  the  region  of  the  s  plane  where  G2  is  the  dominant 
solution,  G2(s,o)  will  have  poles  at  those  points  where  the 
difference  of  the  exponents  c-a-b  is  a  positive  integer. 

b.  Poles  may  arise  in  the  analytic  continuation  of  the  function 
K ( <f> (m)  ,  G2,  ol»  o)  . 

c.  The  function  k(s)  may  have  poles. 

To  study  the  poles  which  arise  according  to  reason  b,  we 
develop  the  integrand  in  powers  of  (oL  -  a) .  In  general  these  will 
be  fractional  powers,  which  depend  upon  s.  Consider  an  expression 
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/ 

r=  rr 


<-l>p(s,«  -  (a"3L)Pts)+1 


0=0, 


(119) 


Evaluating  the  expression  in  the  region  where  the  integral  exists, 
that  is  for  p(s)  =  -1,  one  obtains 

a 

f  <3-oL)p(s)d3  =  (0-aL)p(s)+1  1120) 

0=0 T 

Xj 

The  right  hand  side  can  be  continued  analytically  throughout  the 
s  plane.  A  pole  arises  at  p(s)  =  -1.  Assume  that  this  happens 
at  a  point  s  =  s  .  Then  one  has  in  the  vicinity  of  such  a  pole 

lr 

o 

f  (3-ol>P(S)«  -  jhj  ■  ITpIsWl  +  0<1)  ,121) 

J  p  s=s 

•w  P 

o=oL 

We  have  found  that  the  powers  which  arise  in  K(<p(m),  in  ,  o2,  o) 
are  (oL~o) (5s/3) ~ (7/6)  for  a  =  0,  and  (oL-o) (7s/5) - (3/5)  for 
a  =  1.  Poles  may  therefore  arise  if,  for  a  =  0 


-(7/6)  +  ( 5s/3)  =  -k3,  k3  =  1,  2,  3,... 


that  is  for 


s  =  -(3k3/5)  +  (7/10) 

For  these  values  one  has 

c-a-b  =  k3~l 

The  same  result  is  found  for  a  =  1.  Actually, # the  value  k3  — 
must  be  excluded.  For  then  c-a-b  =0.  In  this  case  §2  =  §3, 
according  to  Eq.  (101)  .  But  we  have  found  that  no  singularity 
will  occur  if  u>  is  formed  with  §3*  In  sxamraary,  poles 
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will  arise  in  the  analytic  continuation  of  K  for  c-a-b  =  1^; 

In  the  above  discussion  we  have  excluded  the  possibility 
that  the  development  of  <j>,  4>t  or  ^  contains  nonintegral  powers 
of  (o-aL) .  If  this  possibility  is  admitted,  then  further  poles 
will  arise  in  the  function  K.  The  particular  solution  for  the 
values  s  pertaining  to  these  poles,  give  individually  the 
propagation  of  a  singularity  along  the  limiting  characteristic. 

This  possibility  is  included  in  the  discussions  of  Reference  4. 

If  <(>,  <J>t  and  4>o  are  regular  at  a=aL,  then  poles  in  the 
function  §9  (reason  a)  and  poles  in  K  occur  at  the  same  values  of  s. 
The  singularities  in  §2  have  the  form  (s-sp  const  §3.  We 
have  seen  above  that  the  functions  §3  do  not  cause  singularities 
in  K.  It  follows  that  no  poles  of  second  order  will  arise  by  the 
coincidence  of  poles  in  §2  and  in  K. 

The  function  k(s)  (Eq.  96))  vanishes  because  of  the  fac¬ 
tor  ( ni-(c-a-b)))  1  for  the  values  of  s  for  which  c-a-b  =  k^.k^  = 
1,2,...,  unless  simultaneously  a  =  ~k2  or  b  =  -k2 ,  k2  =  0,1,2... 

This  factor  cancels  the  effects  of  the  poles  which  occur  in  §2  and 
in  K  ( <}> ,  G2,  aL,  a,  s)  .  If  a  =  0,  then  the  values  of  s,  for  which 
c-a-b  =  k,  and  b  =  -k2  or  a  =  -k2 ,  are  always  different.  The  only 
poles  which  can  arise  are  due  to  k(s) .  It  has  been  demonstrated 
in  the  discussion  following  Eq.  (83),  that  for  a  =  1,  there  are 
some  values  of  s  for  which  simultaneous  c-a-b  =  k^  and  b  =  k2« 

For  these  poles  separate  formulae  must  be  derived. 

Poles  caused  by  the  factor  k(s)  in  Eq.  (114)  arise,  according 

to  Eq.  (96) ,  if  a  =  -k2  or  b  =  -k2,  k2  =  0,  1,  2...  provided  that 

c-a-b  ^  k^,  k^  =  1,  2,  3...  One  recognizes  that  one  has  for  all 

possible  poles  either  a  =  -k2  or  b  =  -k2?  this  implies  that  is 

a  polynomial  in  a.  This  guarantees  that  is  regular  at  a=aL. 

That  particular  solutions  of  only  this  kind  will  be  encountered  must 

be  expected  because  of  the  assumption  that  <j>,  <j>  and  $  are  free 

x  y 

of  singularities  at  a=aT  . 

JLj 
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One  has,  according  to  Eq.  (100) ,  for  these  cases 


^1  _  r(c) r(c-a-b) 
G2  nc-a)r(c-b) 


k^  /  0 ,  1 ,  2 .  .  . 


(122) 


fl  _  ,  U£l 

G2  v  ’  r(a) 


(123) 


These  ratios  can  be  evaluated  for  the  pertinent  values  of  s  =  s^. 
The  specific  values  for  a  =  0  are  found  in  Eqs.  (79)  ,  (80)  and 
(SI).  For  a  =  1,  they  are  computed  from  Eq.  (82),  In  these  cases, 
Eq.  (115)  simplifies  to  the  following  expression 


ip  ^  (s,o)  = 


(m) 


(m) 


(s-sp) (d(k~1)/ds)  G2 


G2  '  (c,s  ) K ( <J) , u>2  f  (t,o,e,sp)  ,0 ,oL 


(124) 


In  the  expression  K,  defined  in  Eq.  117)  ,  we  have  written  u, 
instead  of  the  argument  G2  ; ,  which  might  have  been  used  just 
as  well,  to  emphasize  the  significance  of  the  particular  solutions 
<ii2m^  .  In  these  cases,  the  two  integrals  in  Eq.  (114)  can  be 
combined.  (The  limits  of  the  integral  become  then  0  and  oL .  In 
the  vicinity  of  a  point  s  =  sp  the  o-dependence  of  (s,a)  is 

solely  introduced  by  G2^m\cr,sp).  The  ratio  G^/G2  is  computed 
from  Eqs.  (122)  and  (123) . 

In  cases  where  at  a  point  s  =  Sp  one  has  simultaneously 
b  -  k2,  k2  =  0,  1,  2...  and  c-a-b  =  e  =  k^ ,  k^=l,  2,  3...  one 
obtains  a  finite  limiting  value  for  k(s),  which  can  be  evaluated 
by  Eq.  (98)  .  We  distinguish  between  cases,  where  k2  <  and 
k2  >  kx. 

If  k2  <  k^ ,  then,  according  to  Eq.  (102),  the  function  §2 
remains  finite  in  the  limit  s  ->■  sp.  In  this  limit  §2  consists  of 
the  polynomial  g 2 ,  and  the  function  multiplied  by  a  constant. 
Poles  arise  in  K,  because  of  g2. 
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The  behavior  of  ip  (s,a)  in  the  vicinity  of  sp  is  computed 
by  applying  Eq.  (122)  to  the  expression  (118)  ,  substituting  for 

the  expression  (117)  formed  with  g2>  The  procedure  is  rather 
complicated  because  it  is  based  on  a  series  development  of  the 
integrand.  If  it  should  be  needed,  one  would  probably  carry  it  out 
in  the  physical  plane,  on  the  basis  of  Eq.  (40),  (41)  or  (42).  At 
the  moment  we  simply  write 


K ( o , s )  ~ — “ —  Re sidual  of  K(c^,  oT  ,..»)  at  s  s 

S“Sp  ^  1j  * 

(The  second  limit  of  the  integral  is  unimportant.)  Then,  one  has, 
using  Eq.  (105) 


,  r(c)r(k.) 

'P  (s'a)  =  I^i-  k(sp)  r(k1-k2)rTk1+a)  g21o,sp) 


-(m) 


x  Residual  of  K(<(),  ^  ,  oL,...)  at  s  =  sp ,  k2  <  k^  (125) 


The  function  §2  consists  for  k2  >  kx  according  to  Eq.  (112)  ,  (a)  of 
a  polynomial  ^2  which  contains  powers  of  (1  -  (a/oL))  up  to  k^-1 
with  finite  coefficient ,  (b)  of  a  pole  whose  o  dependence  is  given 

by  the  polynomial  (in  which  k_  is  the  maximum  power  of  (l-o/o^)  , 
and  (c)  an  infinite  series  starting  with  (l-o/aL)  .  A  pole  arises 
because  of  the  contributions  a  and  b.  Invoking  Eq.  (107)  one 
obtains 


♦ 


(m) 


(s,a) 


1 

S-Sp 


k ( Sp) (-)kl~k2 


r(c)ru+k2) 
r (a) r(i+k1)  G3(a'V 


x 


|Residual  of  K(4>, 


=  (m) 
«2 


/ 


o^#»»)  at  s  Sp ) 


1 

de/ds 


r (a+k^) r (l+k2) 
r(a)r(k1)r(i+k2-k1)T 


K(4>, 


(126) 
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SECTION  VII 

THE  SOLUTION  IN  THE  PHYSICAL  PLANE 


In  Eq.  (52)  a  Fourier  decomposition  with  respect  to  0  of  the 
solution  of  the  boundary  value  problem  formulated  in  Section  III 
has  been  made.  The  solution  for  one  Fourier  component,  expressed 
by  means  of  its  Laplace  transform  is  given  by 

+i» 

<Ji*m^(t,0)  =  jjk-  J  (s,o)  exp(st)ds  (127) 

-i°° 


To  ensure  the  behavior  postulated  in  Eqs.  (47)  and  (48),  the  path 
of  integration  must  be  chosen  somewhat  to  the  right  of  the  lines 
Re  s  =  (2/5)  for  a  =  0,  and  Re(s)  =  -(2/7)  for  a  =  1 .  How  far 
this  path  can  be  shifted  to  the  right  depends  upon  the  location  of 
the  neighboring  poles  in  the  s-plane. 

The  integrand  of  Eq.  (127)  depends  upon  s.  The  integrand 

becomes  small  for  negative  t  if  Re  s  is  large  and  positive.  We 

close  the  path  of  integration  by  a  sequence  of  circles  in  the 

right  half  of  the  complex  plane  whose  radius  goes  to  infinity.  It 

is  shown  in  Ref.  4  that  in  a  certain  well  defined  region  in  the 

left  half  of  the  t,a  plane  the  contribution  to  the  integral  due 

to  the  integration  along  these  circles  vanishes  in  the  limit  of 

an  infinite  radius.  In  this  part  of  the  t,a  plane  the  functions 

<f>^m^(t,a)  can  then  be  represented  in  terms  of  the  residuals  at 
/  _  \ 

the  poles  of  '  '  (s,a)  .  With  respect  to  these  poles  one  travels 
in  the  counterclockwise  direction. 

If  the  far  field  conditions  are  satisfied,  then,  according 
to  the  discussions  of  Section  II  following  Eq.  (48),  <}>  ^  =  0  for 
t  <  t  (o)  .  The  region  in  which  the  representation  of  1  in 
terms  of  the  poles  in  the  right  half  of  the  complex  s  plane  is 
valid,  lies  within  the  region  t  <  tg(o).  In  order  for  to  be 

identically  equal  to  zero  for  t  <  t  (o)  it  is  necessary  that 
residuals  at  all  of  these  poles  vanish.  This  observation  is 
decisive  for  the  formulation  of  the  far  field  conditions. 


The  path  of  integration  can  also  be  closed  by  circles  in 
the  left  half  of  the  s  plane  and  the  contribution  of  these 
circles  vanish  for  a  well  defined  region  in  the  right  half  of  the 
t,o  plane.  In  the  present  context  one  is  not  interested  in  a 
representation  of  4> v  ;(t,o)  in  terms  of  the  residuals,  except  for 
the  contribution  of  those  poles,  which  give  a  modification  of  the 
basic  flow  field. 

Expressions  for  i^m^(s,a)  valid  in  the  vicinity  of  the  poles 
have  been  derived  in  Section  VI.  In  the  vicinity  of  a  pole  with 
s  =  they  have  the  form 


s 


c 


i,m 


(m) 

2 


(a,si) 


where  c.  is  a  constant. 

1  ,m 

Closing  the  path  of  integration  by  circles  in  the  right  half  of 
the  complex  s  plane,  one  obtains 


4>(n,)(t,o)  =  -  Y,  ci,m  ^  (a'si}  exP(sit} 


(m) 


(128) 


and 


m= 


4>  ( t 


,o,e)  =  -  ^  cim  G2m)  (a,si'  exP<Bi»t)  exp(ime)  (129) 


m=-°°  i 


The  summation  must  be  extended  over  all  poles  to  the  right  of  the 
lines 


Res  =  2/5  for  a  =  0 
Res  =  -2/7  for  a  =  1 


(130) 


The  individual  terms  are  particular  solutions  of  the  partial 
differential  equation  (37)  ,  which  is  equivalent  to  Eqs.  (12)  or 
(13)  . 
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The  inner  sum  converges,  according  to  Ref.  2  in  that 
region  of  the  t,o  plane  for  which  the  closing  of  the  path  of 
integration  in  the  complex  s  plane  is  permissible..  This  region 
usually  does  not  extend  to  the  curve  t  =  t  (o)  (which  is  the  map 
of  the  outer  surface  s  of  the  computed  flow  field) .  In  order  to 
express  the  postulate  that  <p  =  0  to  the  left  of  t  <  t  (o)  ,  it 

5 

suffices  if  this  condition  is  imposed  in  a  smaller  region.  One 

obtains  the  condition  c.  „  =  0  for  all  n  and  all  poles  to  the 

i,m 

right  of  the  lines  characterized  above.  This  leads  to 


K  (  <j>  ,  w2(m){t,  a,  8,  S;L)  ,  o,  aL)  *  0  (131) 

where  K  is  defined  in  Eq.  (117)  .  The  notation  used  here  is  the 
same  as  in  Eq.  (124) .  K  is  best  evaluated  in  the  physical  plane, 
by  means  of  Eq.  (40)  ,  (41)  or  (42)  .  The  expressions  oi2  are  defined 
in  Eqs.  (118)  .  Here  the  function  §2,  or  also  the  functions 
which  for  these  poles  are  proportional  to  §2,  must  be  substituted. 
They  are  defined  by  Eqs.  (77) . 

One  needs  for  the  evaluation  in  the  physical  plane 

=  y_1(10/3)  (Y+l)“1/3(o/y2)2/5(a+(4/3))_1a  a  =  Q 

(4/3)y(c/y2) (2-a) (a+4/3)"1  a  =  0 

(132) 

y_1(7/10)  <Y+l)”1/3(a/r2)2/7(a+(l/5))-1a  a  =  1 


3a (x,r) 
3r 

=  (2/5)r(o/r2 

) (l-2a) (a+ (1/5) ) ”1 

a  =  1 

One  has 

a/y2  = 

<y«>5/6  p5/2 

(-x)5/2  (2-a)5/2 

a  =  0 

a/ r2  = 

(y+l)7/«  „7/2 

(-x)7/2  (l-2a) 7/2 

(133) 

0  =  1 

3a (x,y)  _ 
3y 

3a (x,r) 

3x 
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These  expressions  are  used  in  a  vicinity  of  o  -  0. 


3o) 

sym 


y^UO/S)  (Y+l)“1/3(a/y2)  (8P/2)  +  (3/10) 
(a+(4/3))“1(aL-a)<5sP/3)_(7/6) 

|[(sp/2)-(VlO)3  (cL-c)-((5sp/3)-(lA))  aj§  (a,sp) 
+  a (aL~a) d^/do| 


=  y-1  (10/3)  (y+l)"1/3(a/y2)  ( (sp/2)  +  (3/10) ) 

(ct+(4/3)  )“1(a1/2(aL-o)  ( /3)  "  (7/6) 

| [  (s  p/2 )  +  ( 2/5 )  ]  ( oL-o)  -  ( ( 5sp/3 )  -  (1/6 ) )  a]  §+a(  aL~a)  (d§/da) } 


Mtosym)/3y  =  y(a/y2)  (Sp/2)  +  (9/1°^aL-o)  (5«p/3) -(7/6)  (a+(4/3))-l 

{t-(10/3)  ( (sp/2)  - (1/10) )  (aL-a)-(esp/3)-(l/6))  (4/3)  (2-a)  ]§(o,Sp) 
+(4/3) (2-o) (oL-o) 09/3o)}  (134) 


8(u)anti)/9y  =  {o/y2)  (V2)  +  (2/5)  (aL-a)  (5sP/3) " (7/6)  (o+(4/3))~1 


|[-  (10/3)  ( (sp/2)-(l/10)  )  (aL-o)o-((5sp/3)-(l/6))  (4/3)  (2 -a)  a 
-(2/3) (2-o) (oL-o) ]§+(4/3) (2 -a) (aL~o) a (d$/dc| 


(3w(m)/3x)  -  y_1(7/10)  {y+l)~1/3om/2(o/z2)  (sP/2i  +  (4/7) 
(a+(l/5))"1(oL-a)  (7sp/5)"(3/5) 

{  [((m/2)  +  (sp/2)  +  (2/7) )  (aL~a)-  ( (7Sp/5)  +  ( 2/5) )  cr]§  (a,  sp) 
+  a(aL-a)  (d^/da)}  exp(-lm0)  (135) 


One  has 

<  \  «  (in)  (m) 

3  a)  m  / 3y  =  —  cos  0  +  im  —  sin  6 

dr  17 

lm\  ■x  (m)  (m) 

3ur  V3z  =  —5 -  sin  0  -  im  - -  cos  0 

<3r  r 


(136) 


Here 

3(jj(m) /3r  =  exP(-im0)a(m/2)-1/2(a/r2)  V2)  +  {11/14) 

(aL-a) (7sp/5)"(3/5) (0+(1/5))-l  (137) 

{[(aL-o) (m/5) (l-2a) - (14/5) ( (sp/2) + (2/7) ) 0)  - 

((7sp/5)  +  (2/5) )  (l-2o)o] §  +  (l-2a)a(aL-CT)ad§/doJ 
for  m  =  0,  this  reduces  to 

^  -  c^o/r2)  (V2>*(11/14)  (0l.0)  (7sp/5)-(3/5)  (o+(1/5))-1 
{[(oL-o)  (-(14/5)  (sp/2 )  +  { 2/7)  )  -  (7Sp/5)  +  (2/5)  (l-2a)]§(0) 

+  (l-2o) (aL-o)d§(0)/dc|  (138) 


s 
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The  values  of  sP  are  given  by  Eq.  (79)  and  (81)  .  For  symmetric 
solutions,  one  has 

sSp  =  (2/5)  +  (6k/5),  k  =  1,  2,...  (139) 

The  value  k  =  0,  which  is  included  in  Eq.  (79)  must  be  omitted 
because  it  lies  outside  of  the  region  given  by  Eq.  (130) ,  it 
amounts  to  a  change  of  the  intensity  of  the  basic  singularity,  and 
such  a  change  must,  of  course,  be  admitted.  For  antisymmetric 
solutions,  one  has 

sp=  1  +  (6k/5)  k  ■  0,  1,  2,...  (140) 

The  values  of  s  for  a  =  1  are  computed  from  Eq.  (82)  .  The 
value  for  m  =  0  and  k  =  0  must  be  omitted  because  they  lie  outside 
(although  at  the  border)  of  the  region  characterized  by  Eq.  (130) . 
According  to  Eq.  (100) ,  which  is  valid  for  these  poles,  the 
function  §2  is  proportional  to  As  long  as  the  pole  strength 

is  set  equal  to  zero,  the  factor  of  proportionality  is,  of  course, 
unimportant.  It  is  best  to  compute  §2  from  the  first  of  Eqs.  (77) . 

There  exists  the  possibility  to  extend  the  present  approach 
to  boundary  conditions  for  a  free  stream  Mach  number  close  to  one. 

In  such  cases  one  would  need  the  actual  strength  of  the  poles. 

This  (and  other  factors)  have  therefore  been  carried  in  the  analysis. 

Assume  that  a  correction  to  the  flow  field  has  been  computed 
with  the  use  of  the  far  field  conditions  Eqs.  (131) .  As  a  second 
step,  one  adjusts  (if  necessary)  the  intensity  and  the  origin  for 
the  basic  singularity.  This  is  done  on  the  basis  of  Eqs.  (89) 
and  (90)  ,  which  express  the  effect  of  a  change  of  y  and  of  the 
location  of  the  origin  by  means  of  certain  particular  solutions. 

For  this  purpose,  one  closes  the  curve  of  the  original 
inversion  integral  by  circles  in  the  left  half  of  the  complex  s 
plane.  One  then  obtains 

*{m)(t,0)  =  f  cim  G2m)  (cT'si)  exP<sit^  (141) 
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The  summation  is  over  all  poles  that  lie  to  the  left  of  the  line 
given  by  Egs.  (130)  or  on  this  line.  The  expression  is  slightly 
different  for  poles  where  simultaneously  c-a-b  =  ,  k1=l,...  and 

b  =  -k2f  k2  =  0,  1,  2,...  But  in  the  present  discussion  these 
poles  are  of  no  interest.  The  values  of  c.  „ equal  to  the  residuals 
at  the  individual  points,  which  are  given  by  Eq.  (124) ,  (125)  and 
(126)  .  The  adjustment  of  the  basic  solution  is  based  on  Eqs.  (89) 
or  (90) . 

This  requires  for  a  *  0  that  we  find  the  contributions  of 

Gsym  (a'2/5)'  Gsym^a  '~2/5)  and  Ganti  •  Tlle  corresponding 

contributions  for  a  =  1  are  G (a,  -2/7),  g(0)(o,-( 6/7)),  and 

G^)(o,  -9/7).  Let  us  denote  the  constants  c.  „  in  the  same 

i  ,m 

manner  as  the  corresponding  functions  G.  These  constants  can  be 
evaluated  numerically.  Combining  these  expressions  with  the 
basic  flow,  one  obtains  an  expression  (for  a  =  0) 

i0  -  y2/5  r<C>  +  csym,  2/5  Gsym(o'  2/5>  exP(2t'5> 

+  csym, -2/5  G3ymll,'-2/5)  exPl-2t'5> 

+  cantt,-3/5  Ganti "3/5)  e*P<-3t/5> 

According  to  Eq.  (89)  this  is  an  approximation  to 

$0  =  (y  -  Ay) 2/5  (y  +  Ay)3  f(f)  (142) 

fc/ith  r  defined  in  Eq.(85))  if  one  sets 

ip  =  (1/5) y-2  csym2/| 

4x  -  (Y  +  1)1/3  p"2  csymr_2/s  (143) 

4y  -  -2p'3  cantl,.3,5 

A  similar  procedure  is  applied  for  a  *  1. 
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In  the  following  iteration  step  the  expression  <j>  is  then 
used  for  the  basic  flow  field.  This  correction  needs  to  be  carried 
out,  only  if  Ay,  Ax,  and  Ay  are  fairly  large.  Sometimes  it  may 
suffice,  if  one  makes  a  correction  only  for  the  dominant  term,  which 
is  given  by  Ay. 


SECTION  VIII 

NUMERICAL  IMPLEMENTATION 


It  is  assumed  that  an  approximation  to  the  flow  field  is 
available  and  that  one  wants  to  improve  it  iteratively.  One 
needs  for  the  application  of  the  far  field  a  characterization  of 
the  basic  field  for  the  outer  region  by  the  value  of  y  and  the 
coordinates  of  the  origin  to  which  it  refers.  The  correction  to 
the  flow  field  is  then  computed  using  Egs.  (131)  as  conditions 
to  be  imposed  at  the  distant  boundary  of  the  computed  flow  field. 

For  the  field  so  obtained  determines,  in  a  second  step,  corrections 
to  the  basic  field. 

The  computed  part  of  the  flow  field  will  frequently  consist 

of  one  portion  which  is  given  by  the  result  of  the  previous 

iteration  and  a  correction  to  it,  which  is  to  be  determined.  Let 

$ .  _  be  the  sum  of  these  two  contributions.  The  outer  field 

inner 

is  given  according  to  Eq.  (11)  as  the  sum  of  the  basic  field 
and  a  correction  <}>.  $>Q  is  defined  by  Eqs.  (6)  ,  (7)  ,  and  (10)  . 

The  function  f(^)  is  given  in  parametric  form  in  Eqs.  (27)  and  (28) 
for  a  =  0  and  a  =  1,  respectively.  The  position  of  the  origin  with 
respect  to  which  c(via  x,y(and  z) )  is  defined  may  vary  from 
iteration  to  iteration.  The  conditions  (131)  refer  to 

♦  =  ^ inner  '  *0  (144) 

In  the  application  of  Eqs.  (.131)  one  must  therefore  have  a 

characterization  of  4>. _ _ (.and  itB  gradient),  for  instance  by  the 

inner 

values  of  <j)  at  a  system  of  grid  points,  or  in  a  finite  element 

characterization.  Besides  one  needs  the  function  and  the 

(m)  y 

"weight  functions" 

Assume  that  the  outer  boundary  S  of  the  computed  flow  field 
is  defined  by  the  coordinates  at  a  finite  number  of  points.  The 
functions  oT  have  the  form  G(C,s)  exp  (st)  =  G(C,s)p8.  The 
function  G(C,s)  oscillates  in  the  subsonic  region  with  an  amplitude 
which  does  not  change  much  even  if  s^®.  If  along  the  outer  boundary 


of  the  computed  region  p  changes  considerably,  then  there  is  a 
danger  that  the  functions  w  for  large  values  of  8  are  close  to 
linear  dependence.  (For  large  values  of  s  the  functions  w  are 
then  important  only  at  those  points  where  p  is  large.)  It  is 
therefore  desirable  to  choose  the  outer  surface  of  the  computed 
flow  field  close  to  a  surface  p  =  const  (see  Figure  1) . 

The  points  (x,y) ,  (x,r)  or  (x,y,z)  used  in  the  computation 
of  the  inner  field  are  considered  as  fixed.  A  subscript  0 
characterizes  the  coordinates  of  the  current  origin  used  for  the 
representation  of  the  outer  flow  field. 

One  computes  from  Eq.  (7)  for  the  chosen  points  of  S 


C  =  (y+1)  **1//3  (x-x0)  (y~yQ)  ~4/,S  for  a  =  0 


where 


C  =  (y+l)-1/3(x-x0)  iT4/7 


f  -icy  -  yn) 2  +  (z  -  zn)2]i/2 


(145) 


for  a  =  1 


or  in  axisyiranetric  flow. 


One  has,  furthermore 


Next,  one  determines 


r  =  r 


(146) 


=  arctg 


y_1c 


and  then  finds  a,  either  from  a  table,  or  by  a  solution  of  the 
second  of  Eqs.  (27)  or  (28)  ,  for  instance  by  means  of  a  Newton 
iteration.  This  then  allows  one  to  determine  t  from  the  second  of 

3  _ 

Eqs.  (27)  or  (28)  and  f  ■  y  f  from  Eq.  (10) .  Then  one  obtains 
from  Eq.  (6) . 
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Figure  1.  The  £-p  System  for  Axisymmetric  Flows. 


I 


$0  =  (y  -  yQ) 2/5  f  for  a  =  0 
$q  =  r  f  for  a  =  1 


These  expressions  are  used  in  Eq.  (144)  . 

The  integrals  (131)  will  probably  be  evaluated  in  the  form 
of  Eqs.  (40),  (41),  and  (42)  for,  respectively,  three-dimensional 
and  plane  flow  in  Cartesian  coordinates,  and  three  dimensional 
flow  in  cylindrical  coordinates.  Axisymmetric  flows  arise  by 
specialization  of  Eq.  (42) .  The  coordinates  used  in  these 
systems  are  those  used  for  the  computation  of  the  inner  flow  field. 
The  functions  to  to  be  used  are  computed  from  Eq.  (116)  in  the 
coordinates  of  the  outer  field.  The  functions  are  defined  in 
Eqs.  (118) .  The  function  §  occurring  in  Eqs.  (118)  are  defined 
in  Eq.  (77)  in  terms  of  a.  Here  oT  =  16/3  for  a  =  0,  and  ar  =  6/5 
for  a  =  1.  The  values  of  a,  b,  and  c  occurring  in  these  formulae 
are  found  in  Eq.  (76)  ,  the  value  of  R,  which  occurs  here  is  given 
in  Eq.  (72) .  The  values  of  s  for  which  these  functions  ^  are  to 
be  given  for  the  plane  case  in  Eqs.  (139)  and  (140)  .  For  the 
three-dimensional  case,  they  are  computed  from  Eq.  (82)  ,  with  the 
provision  taken  from  Eq.  (13)  that  Re(s)  >  -  2/7.  The  number  of 
functions  w  ^  to  be  used  equals  the  number  of  parameters  which 
one  uses  to  characterize  the  function  4>  at  the  outer  boundary 
surface  s. 


The  formulation  of  the  conditions  Eq.  (131)  for  a  sufficient 

number  of  functions  gives  as  many  linear  relations  of  a 

global  character  between  the  4> ,  4  ,  and  <p  at  the  outer  boundary 

x  y 

as  are  necessary  for  the  computation  of  the  flow  field. 


After  the  values  of  $,  and  are  found .  one  uses  Eqs*  ..(1431 


and  (142)  to  find  a  readjusted  basic  field  4>0  which  is  to 


be  used  in  the  next  iteration  step.  The  constants  c 


(m) 


sym,2/5' 


(s,a)  at  the  respective  poles. 


c _  0 /c  are  the  residuals  of  \b 

sym  f  "*«/  5 

They  are  found  from  Eq.  (124)  with  the  definitions  Eqs.  (122)  or 


(123) . 
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The  values  of  s  needed  for  this  adjustment  are  (2/5)  and  -2/5) 
for  the  symmetric  solutions  and  (-3/5)  for  antisymmetric  solutions 
in  the  plane  case.  For  the  case  a  =  1,  the  values  of  s  are  taken 
from  Eq.  (90)  ,  they  are  s  =  -2/7  and  s  =  -6/7  for  m  =  0,  and 
s  =  -9/7  for  m  =  1. 

In  evaluating  the  integrals  K,  one  must  observe  the  behavior 
of  the  integrand  at  the  limiting  characteristic  a  =  aL,  which  is 
determined  by  the  behavior  of  a».  One  recognizes  from  Eq.  (118) 
that  for  a  =  aT  and  y  regular 

,  .  (5s/3)-(l/6)  *  - 

a)  ~  (a-  -  a)  '  for  a  =  0 

,  .  (7s/5)  +  ( 2/5)  . 

w~(aT-a)  '  '  /  for  a  =  i 

Li 

One  obtains,  in  particular 


0J  " 

(aL  ‘ 

-  O  1/2 

s  = 

2/5 

0)  ~ 

(°L  ’ 

■  a, ~5/6 

s  = 

-2/5 

j  a  «  0 

0)  ~ 

(aL  ’ 

-  c)-7/6 

s  = 

-3/5 

1 

u  ~  (aL  -  a)0  for  s  =  -2/7 

u  ~  (aL  -  a) _4/5  for  s  =  -6/7  a  =  1 

w  ~  (cL  -  a)  7^5  for  s  =  -9/7 

Notice  that  in  the  integral  K  first  derivatives  of  u  are  encountered. 
Except  for  the  integrands  for  s  =  2/5,  (a  =  0)  for  which  the 
integral  K  is  defined  as  an  improper  integral  and  for  s  =  -2/7  for 
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a  =  1,  for  which  no  singularity  is  encountered,  the  integrals 
are  defined  by  analytic  continuation  (Eq.  (121)). 

One  needs  for  this  purpose  the  first  few  terms  of  a  power 
series  development  of  0  and  its  gradient  in  terms  of  (a  -  a or 
an  equivalent  quantity.  They  must  be  extracted  from  the  numerical 
data  for  0.  Only  a  few  of  these  terms  can  be  obtained  with 
acceptable  accuracy.  The  integral  K  is  evaluated  by  one  of  the 
familiar  procedures  for  0  <  a  <  a1  oL  and  by  means  of  formulae 
developed  from  Eq.  (121)  for  the  region  a^<  a  <  a^. 

For  the  values  of  s  which  appear  in  Eq.  (131)  the  integrals 
are  integrable.  It  may,  however,  be  worthwhile  to  examine  to  what 
extent  the  behavior  of  the  function  w  at  o  =  aL  must  be  taken  into 
account.  These  are  the  main  considerations  needed  for  a  practical 
application. 
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SECTION  IX 
CONCLUDING  REMARKS 

The  far  field  conditions  (Eq.  (131))  are  global  conditions; 
that  is,  they  give  relations  connecting  all  points  of  the  surface  S. 
This  makes  their  application  in  conjunction  with  an  iterative 
procedure  (as  for  instance  the  Murman  Cole  iterations)  rather 
difficult.  They  are  rather  well  suited  for  a  correction  to  an 
existing  field  by  means  of  a  Newton-Raphson  iteration  in  combination 
with  a  direct  elimination  procedure  in  the  resulting  linear  equations. 
In  such  a  procedure  the  sequence  of  the  elimination  of  the  unknown 
is  important.  One  best  arranges  the  problem  in  a  manner  that  one 
obtains  a  large  block  multidiagonal  matrix.  The  vector  of  the 
unknown  is  then  partitioned,  the  first  of  the  subvectors  which 
arise  in  this  manner  express  data  at  the  surface  of  the  body  for 
which  the  flow  is  computed,  the  last  of  the  subvectors  gives  data 
at  the  surface  S.  Usually  the  individual  blocks  are  banded  matrices, 
while  the  matrices  derived  from  Eq.  (131) ,  which  form  the  last  row 
of  blocks,  are  full.  This  does  not  matter  because  in  the 
elimination  process  for  the  subvectors,  one  obtains  full  matrices 
rather  quickly.  These  matrices  have  the  same  size  as  the  matrices 
expressing  Eq.  (131) .  In  any  case,  it  seems  desirable  to  use  a 
coarser  mesh  at  larger  distances  from  the  body. 

The  conditions  Eqs.  (131)  are  exact,  provided  that  the 
linearization  in  the  outer  field  used  in  Eq.  (11)  is  applicable. 

The  global  character  of  these  particular  solutions  is  unavoidable. 
Reference  2  derives  for  another  problem  approximate  conditions, 
valid  for  large  distances,  which  have  a  local  character.  They 
would  then  be  applicable  for  an  iteration  of  the  Murman-Cole  type. 

The  principal  idea  is  to  set  up  relations  which  are  compatible  with 
the  dominant  terms  at  infinity.  The  dominant  term  in  the  basic 
flow  is 

<j>  =  y2//5  f(C)  for  a  =  0 

<p  =  r”2/7f(c)  for  a  -  1 
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Along  a  line  c  =  const  one  has 


d<t>  _  2  1 
3y  5  y 


0 


for  a  =  0 


|i+2i^=0  for  a  -  1 

These  are  localized  for  field  conditions.  The  first  perturbation 
occurs  for  a  =  0,  at  s  =  0  according  to  Eq.  (80).  This  leads  to 
a  boundary  condition 


<f7+7>  'I?  0 

(The  derivatives  are  to  be  formed  along  lines  x,  -  const.) 

Indeed,  if  4>  =  y2/^k^(5),  then  the  condition  is  satisfied  because 
of  the  second  operator;  if  one  substitutes  <p  =  k2(c),  then  the 
second  operator  gives  -2/5 ( 1/y) k2 ( c) ,  which  gives  zero  if  the  first 
operator  is  applied  to  it.  Here  kjJOand  k2(c)  are  arbitrary 
functions  of  £.  The  second  derivative  with  respect  to  y  must  then 
be  expressed  by  derivatives  along  the  surface  S  using  the  partial 
differential  equations  for  *. 

In  a  similar  manner,  one  can  obtain  improved  conditions  for 
the  three-dimensional  case.  Conditions  of  this  kind  are,  of  course, 
rather  attractive,  and  they  are  preferable  to  such  assumptions  as 
<J>  =  a*x  (a  sonic  free  jet)  or  <pw  =  0  at  a  distant  line  which  actually 
means  a  rigid  wind  tunrfel  wall. 

The  use  of  the  rigorous  conditions  would  allow  one  to  reduce 
the  size  of  the  computer  flow  field.  Whether  this  is  worthwhile 


from  the  point  of  view  of  computational  economy  remains  to  be  seen. 

The  present  article  provides  the  background  for  making  this  comparison. 
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